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ABSTRACT 

We analyze basic properties of about 6000 Lyman-a absorbers observed in high resolu- 
tion spectra of 19 quasars. We compare their observed characteristics with predictions 
of our model of formation and evolution of absorbers and dark matter (DM) pancakes 
and voids based on the Zel'dovich theory of gravitational instability This model as- 
serts that absorbers are formed in the course of both linear and non linear adiabatic 
and shock compression of DM and gaseous matter. Our model is consistent with sim- 
ulations of structure formation, describes reasonably well the Large Scale Structure 
(LSS) observed in the distribution of galaxies at small redshifts and emphasizes the 
generic similarity of the process of formation of LSS and absorbers. Using this model 
we are able to link the column density and overdensity of DM and gaseous component 
with the observed column density of neutral hydrogen, redshifts and Doppler param- 
eters of absorbers. We show that the colder absorbers are associated with rapidly 
expanded underdense regions of galactic scale. 

We extend the method of measuring the power spectrum of initial perturbations 
proposed in Demianski & Doroshkevich (2003b). The observed separations between 
absorbers and their DM column density are linked with the correlation function of the 
initial velocity field. Applying this method to our sample of absorbers we recover the 
CDM like power spectrum at scales 10ft- -1 > D > 0.15ft _1 Mpc with a precision of 
~ 15%. However at scales ~ 3 — 150/i _1 kpc the measured and CDM-like spectra are 
different. This result suggests a possible complex inflation with generation of excess 
power at small scales. Both confirmation of the CDM-like shape of the initial power 
spectrum or detection of its distortions at small scales are equally important for the 
widely discussed problems of physics of the early Universe, galaxy formation, and 
reheating of the Universe. 

Key words: cosmology: large-scale structure of the Universe — quasars: absorption: 
general — surveys. 



1 INTRODUCTION 

One of the most perspective methods to study the processes 
responsible for the formation and evolution of the structure 
in the Universe is the analysis of properties of absorbers 
observed in spectra of the farthest quasars. The great po- 
tential of such investigations was discussed already by Oort 
(1981, 1984) just after Sargent et al. (1980) established the 
intergalactic nature of the Lyman-Q forest. Indeed, the ab- 
sorption lines trace the small scale distribution of hydrogen 
along the line of sight at redshifts z > 2 when matter is not 
yet strongly clustered and its observed characteristics can be 



more easily interpreted. The available Keck and VLT high 
resolution observations of the Lyman- a forest provide a rea- 
sonable database and allow one to apply statistical methods 
for their analysis. 

The composition and spatial distribution of the ob- 
served absorbers is complicated. Thus, at large redshifts the 
population of rich metal systems including Ly-damped and 
Ly-limit systems are rare and majority of observed absorbers 
is associated with isolated low mass HI clouds. At low red- 
shifts a significant number of stronger Ly-a lines and metal 
systems is associated with galaxies (Bergeron et al. 1992; 
Lanzetta et al. 1995; Tytler 1995; Le Brune et al. 1996). 
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However as was recently shown by Penton, Stock and Shull 
(2000; 2002) and McLin et al. (2002), even at small red- 
shifts some absorbers are associated with galaxy filaments 
while others are found within galaxy voids. These results 
suggest that the population of weaker absorbers dominating 
at higher redshifts can be associated with weaker structure 
elements formed by the non luminous baryonic and DM com- 
ponents. They also suggest that the Ly-a forest can be con- 
sidered as a low mass component of the generic Large Scale 
Structure (LSS) which is seen in simulated and observed 
spatial matter distribution. This means that absorbers at 
high z trace the DM structure which is qualitatively similar 
to the rescaled one observed at small redshifts. 

In this paper we investigate a sample of ~ 6000 ab- 
sorbers observed in 19 high resolution spectra of QSOs and 
compare their properties with the improved model of ab- 
sorbers proposed in Demianski & Doroshkevich (2003 a,b, 
hereafter Paper I & Paper II). We assume that absorbers are 
dominated by long-lived gravitationally bound and partially 
relaxed clouds composed of both DM and baryonic compo- 
nents. The fact that we can observe galaxies and quasars 
at 2 > 3 demonstrates the existence of strong density per- 
turbations already then. Here we show that at these red- 
shifts there are also strong negative density perturbations 
of a galactic scale which can be identified with rapidly ex- 
panded underdense regions. Such regions are naturally as- 
sociated with the colder absorbers. 

Our model explains the self similar character of evolu- 
tion of the observed Doppler parameter, HI column density 
and absorber separation. Such evolution implies that the 
mean values of these characteristics slowly vary with red- 
shift while their probability distribution functions (PDFs) 
remain unchanged. This model links the observed and other 
physical characteristics of absorbers - such as their DM col- 
umn density, size, and fraction of matter associated with 
absorbers - and allows us also to identify several subpopu- 
lations of absorbers with different evolutionary histories. 

We treat the evolution of structure as a random pro- 
cess of formation and merging of Zel'dovich pancakes, their 
transverse expansion and/or compression and successive 
transformation into filaments and high density clouds. Later 
on the hierarchical merging of pancakes, filaments and 
clouds forms rich galaxy walls observed at small redshifts. 
Impact of these factors is clearly seen in high resolution nu- 
merical simulations of evolution of the LSS (for review see, 
e.g. Frenk 2002). 

Theoretical expectations of our model are based on the 
Zel'dovich approximate theory of gravitational instability 
(Zel'dovich 1970; Shandarin & Zel'dovich 1988). As is well 
known, it correctly describes only the linear and weakly 
nonlinear stages of the structure formation and cannot de- 
scribe later stages of evolution of structure elements. In spite 
of this, the statistical approach proposed in (Demianski & 
Doroshkevich 1999, 2004a; hereafter DD99 & DD04) nicely 
describes the main properties of observed and simulated LSS 
(Demianski et al. 2000; Doroshkevich, Tucker, Allam & Way 
2004) without any smoothing or filtering procedure. 

Presently various observations are used to determine 
the power spectrum of the initial density perturbations. Its 
amplitude and its shape on scales > 10ft _1 Mpc are approx- 
imately established by investigations of the microwave relic 
radiation (Spergel et al. 2003, 2006) and the structure of the 



Universe at z <C 1 detected in large redshift surveys (Perci- 
val et al. 2001; Tegmark, Hamilton & Xu 2003; Verde et al. 
2002, 2003) and weak lensing data (see, e.g., Hoekstra, Yee 
& Gladders 2002). 

The shape of the initial power spectrum on scales lO/i -1 
Mpc - l/i _1 Mpc can be tested at high redshifts where it is 
not yet strongly distorted by nonlinear evolution (Croft et 
al. 1998, 2002; Nusser & Haehnelt 2000; Gnedin & Hamilton 
2002; Viel et al. 2004 a,b; Kim et al. 2004; McDonald et al. 
2004; Seljak et al. 2004; Zaroubi et al. 2005). The method 
used in these papers is surprisingly universal and is success- 
fully applied to spectra observed with both high and mod- 
erate resolution. It utilizes the measured transmitted flux 
only and does not require preliminary determination of col- 
umn density, Doppler parameters and even discrimination 
of hydrogen and metal line systems (McDonald et al. 2004, 
Seljak et al. 2004). In spite of this, it successfully restores 
the CDM-like power spectrum down to scales ~ Ih^ 1 Mpc. 
Recent results on reconstruction of the initial power spec- 
trum are summarized and discussed in many papers (see, 
e.g., Tegmark and Zaldarriaga 2003; Wang et al. 2003; Zal- 
darriaga, Scoccimorro & Hui 2003; Peiris et al. 2003; Spergel 
et al. 2003,2006; Tegmark et al. 2004; McDonald et al. 2004, 
Seljak et al. 2004). 

A straightforward method of reconstruction of the ini- 
tial power spectrum from the observed characteristics of ab- 
sorbers was proposed and tested in Paper II. This method 
can be used to recover the initial power spectrum down to 
unprecedentedly small scale. In contrast with previous inves- 
tigations (Croft et al. 1998, 2002; Nusser & Haehnelt 2000; 
Viel et al. 2004b; McDonald et al. 2004) we analyze the sepa- 
ration between adjacent absorbers and their column density 
rather than the flux or smoothed density field. This means 
that our results are not restricted by the standard factors 
such as the Nyquist limit, the impact of nonlinear processes, 
the unknown matter distribution between absorbers or their 
peculiar velocities. This approach successfully complements 
investigations of the power spectrum mentioned above. 

Here we improve the analysis discussed in Paper II by 
using a richer observed sample and a more refined model of 
absorbers. We use two independent methods of determina- 
tion of the initial power spectrum. The first one is based on 
measurements of the separation between adjacent absorbers, 
while the second one, proposed in Paper II, uses measure- 
ments of the column density of absorbers. Both approaches 
allow one to determine the spectrum down to the scale of 
~ 5 - 10ft _1 kpc. At scales (0.15 - 10)/i _1 Mpc our results 
coincide with those expected for the CDM-like power spec- 
trum and Gaussian perturbations with the precision of ~ 
15%. However, we have found some evidence that at scales 
< 0.15ft _1 Mpc the initial power spectrum differs from the 
CDM-like one suggesting a complex inflation with genera- 
tion of excess power at small scales. Such excess power accel- 
erates the process of galaxy formation at high redshifts and 
can shift the epoch of reheating of the Universe to higher 
redshifts. 

At present we have only limited information about the 
properties of the background gas and the UV radiation 
(Haardt & Madau 1996; Scott et al. 2000, 2002; Schaye 
et al. 2000; McDonald & Miralda-Escude 2001; McDonald 
et al. 2000, 2001; Theuns et al. 2002 a, b; Levshakov at 
al. 2003; Boksenberg, Sargent & Rauch 2003; Demianski & 
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Doroshkevich 2004b) and therefor some numerical factors 
in our model remain undetermined. This means that our 
approach should be tested on representative numerical sim- 
ulations that more accurately follow the process of forma- 
tion and disruption of pancakes and filaments and provide 
a unified picture of the process of absorbers formation and 
evolution (see, e.g., Weinberg et al. 1998; Zhang et al. 1998; 
Dave et al. 1999; Theuns et al. 1999, 2000). But so far such 
simulations are performed mainly in small boxes what re- 
stricts their representativity, introduces artificial cutoffs in 
the power spectrum and complicates the quantitative de- 
scription of structure evolution (see more detailed discussion 
in Gnedin & Hamilton 2002; Tcgmark & Zaldarriaga 2002; 
Zaldarriaga Scoccimoro & Hui 2002; Seljak, McDonald & 
Makarov 2003; Manning 2003; Paper II). As was shown by 
Meiksin, Bryan & Machacek (2001), the available simula- 
tions reproduce quite well the characteristics of the flux but 
cannot restore other observed characteristics of the forest. 
This means that first of all numerical simulations should be 
improved (see more detailed discussion in Paper II). 

Comparison of results obtained in Paper I and Paper 
II and in this paper demonstrates that the quality and rep- 
resentativity of the sample of observed absorbers are very 
important for the reconstruction of processes of absorbers 
formation and evolution. Thus, richer sample of the observed 
absorbers makes it possible to select and investigate several 
representative subsamples of absorbers with different evo- 
lutionary histories. None the less, our analysis indicates a 
possible deficit of weaker absorbers and pairs of close ab- 
sorbers what in turn could be related to insufficient sensi- 
tivity of the process of absorbers' identification. Thus, the 
number and parameters of absorbers identified for the same 
quasar depend upon the used identification procedure. Fur- 
ther progress can be achieved first of all with richer samples 
covering the range of redshifts at least up to z ~ 5 what 
would allow one to perform a complex investigation of the 
early period of structure evolution. 

This paper is organized as follows. In Sec. 2 the obser- 
vational databases used in our analysis are presented and 
statistical characteristics of the observed parameters of ab- 
sorbers are obtained. The theoretical models of the structure 
evolution are discussed in Sec. 3. The results of statistical 
analysis of the model dependent parameters and the derived 
initial power spectrum are given in Sec. 4. Discussion and 
conclusions can be found in Sec. 5. 



2 OBSERVED CHARACTERISTICS OF 
ABSORBERS 

2.1 Properties of the homogeneously distributed 
matter 

In this paper we consider the spatially flat ACDM model of 
the Universe with the Hubble parameter and mean density 
given by: 

h 2 ( z ) = ff 2 «m(i + zf[i + n A /n m (i + z)- a ] , 



(n b (z)) = 2.4 • 10~ 7 (1 + z) 3 (Q b h 2 /0.02)cm~ 3 , 
{p m {z)) = ^Sl m (l + z) 3 , H = 100ft km/s/Mpc. 



(1) 



Here fl m = 0.3&£7a = 0.7 are the dimensionless matter 
density and the cosmological constant (dark energy), Q b is 
the dimensionless mean density of baryons, and h = 0.7 is 
the dimensionless Hubble constant. 

Properties of the compressed gas can be suitably related 
to the parameters of homogeneously distributed gas, which 
were discussed in many papers (see, e.g., Ikeuchi & Ostriker 
1986; Haardt & Madau 1996; Hui & Gnedin 1997; Scott et 
al. 2000; McDonald et al. 2001; Theuns et al. 2002a, b; 
Demiahski & Doroshkevich 2004b). In this paper we con- 
sider evolution of absorbers at observed redshits z < 4 when 
weak variations of the gas entropy are determined by the 
interaction of the gas with the UV background. Thus, the 
expected intensity of the UV background radiation, can be 
fitted as follows: 



J(z, v) = J21 



10" 



erg 



s cm 2 sr Hz ' 



(2) 



where vh = 3.3 • 10 Hz, a 7 w 1.5 and the dimensionless 
factor J2i(z) describes redshift variations of the intensity. 
The mean temperature, {T bg }, of homogeneous gas can be 
taken as 



(T bs ) 



3.5* 4 6/7 etf-io 4 ^, 



(3) 



66 



fl b h 2 3.5 / 0.15 \ 3/4 



0.02 2 + a 7 \Vt m h 2 



24 = 



where fcs & mu are the Boltzmann constant and the mass of 
the hydrogen atom and a 7 is the power index of spectrum 
of the ionizing background in (2). Analyzing the observed 
characteristics of absorbers McDonald et al. (2001) estimate 
the background temperature as T bg w (2 ± 0.2) • 10 4 if at 
z ~ 2, what is close to (3). 

At this period the gas entropy can be characterized by 
the function 



(T bg ) 



W^ 2/3 keV.cm 2 . (4) 



n b h 2 



As was shown in Demiahski & Doroshkevich (2004b), the 
function F bg oc {n b / {n b )) 0A only weakly depends upon varia- 
tions of the expansion rate. This means that the background 
temperature and Doppler parameter vary as 

T bg oc (n b / (n b )) 2 ' 3 , b bg oc (n b j (n b }) 1/3 , (5) 

and the mildly nonlinear compression or expansion of mater 
occurs almost adiabatically. Thus, within compressed or 
slowly expanded regions with n b > (n b ) we can expect 
that T bg > (T bg ) while in rapidly expanded regions with 
n b <§C (n b ) we can expect that T bg <JC (T bg ) . 

Under the assumption of ionization equilibrium of the 

gas, 

— ) —' (6) 



(n H ) _ a r {n b ) 
(n b ) r 7 



4.4- 10" 



where a r (T) is the recombination coefficient (Black, 1981) 
and T 7 characterizes the rate of ionization by the UV back- 
ground, the fraction of neutral hydrogen is 

x bg = {n H )/{n b } = £ (1 + z) 3i/14 , (7) 
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Table 1. QSO spectra used in our analysis 
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2002 4. Hu ct al., (1995), 5. Djorgovski ct al. (2001) 6. Kirkman 
& Tytler (1997), 7. Cristiani & D'Odorico (2000), 
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1.2 ■ 1(T 7 £Vi 2 



, r 7 (^) = io- 12 ri 2 («)s- 1 , 



T 12 (z) 0.02 69 
ri 2 (z) = 12.6J 2 i(«)(3 + a 7 )- 1 . 

The redshift variations of the rate of ionization, Ti2(z), 
produced by the UV radiation of quasars were discussed by 
Haardt & Madau (1996) and later on tested and corrected 
by Demianski & Doroshkevich (2004b) with large observed 
sample of QSOs. For a-, w 1.5 the expected ionization rate 
of hydrogen is fitted by the expression: 



ri2~7cxp[-(z-2.35) 2 /2] 



(8) 



Variations of the power index a 7 with time and space gen- 
erate random variations of Ti2. With these parameters of 
background we have for the Gunn-Peterson optical depth 



TGP(Z) 



0.34 



r 12 e 



3/7 
bg 



27/7 



0.02 



0.15 



(9) 



tgp(2) « 0.015, r GP (5)^8.2. 



However these estimates of Ti 2 should be corrected for 
absorption and reemission of radiation within high density 
clouds (Haardt & Madau 1996). This effect decreases the 
ionization rate of both hydrogen and helium and enhances 
a possible spatial variations of intensity of the UV back- 
ground. According to the observational estimates of Scott et 
al. (2002) Ti2 ~ 1-4 at z w 2 (see also Levshakov et al. 2003; 
Boksenberg, Sargent & Rauch 2003). On the other hand, Eq. 
(8) takes into account only the contribution of QSOs and, 
so, it underestimates the ionization rate at both z < 2 and 
z > 3 where contribution of other probable sources of radi- 
ation becomes important. 



2.2 The database. 

The present analysis is based on 19 high resolution spectra 
listed in Table 1. These spectra contain 7 770 absorbers. For 
further discussion we selected the sample of 7 430 absorbers 
with 11.9 < IgNni < 15 and 7411 distances between neigh- 
boring absorbers. This sample will be used for discussion 
of the correlation function of initial velocity field. For de- 
tailed investigation we selected a more homogeneous sample 
of 6 270 absorbers and 6 251 separations with b > 5km/s 
and we restricted errors of measurement by the conditions 
A lg Nhi < 0.2 and A6 < 0.36. The chosen low limit of b is 
close to the spectral resolution in this sample. To test the 
sample dependence of the correlation function of initial ve- 
locity field we use for comparison the sample of 14 QSOs 
with 4 036 absorbers investigated in Paper I and Paper II . 

However, a list of absorbers depends also upon the 
method of line identification and, for example, for QSO 
0636+680 two spectra listed in Table 1 include different 
number of absorbers. This example shows that the meth- 
ods of line identification should be unified and improved. 
None the less, dispersions of absorbers characteristics dis- 
cussed below are defined mainly by their broad distribution 
functions and by the completeness of the samples. Because 
of this, in this paper we discuss the scatter of only the more 
interesting quantitative characteristics of absorbers. 

As is seen from Fig. 1, the redshift distribution of ab- 
sorbers is non homogeneous and the majority of absorbers 
are concentrated at 2 < z < 3.5 . This means that some of 
the discussed here characteristics of absorbers are derived 
mainly from this range of redshifts. Absorbers at z > 3.5 
were identified mainly in spectra of QSOs 0000-260 and 
1055+461. In this range and at z < 2 the statistics of lines 
is not sufficient. 



2.3 Observed characteristics of absorbers 

For the sample of 6 270 absorbers the redshift variations 
of the three mean observed characteristics of absorbers, 
namely, the Doppler parameter, (6), the column density of 
neutral hydrogen, (IgN^j) = (\g{NHi / z\)) , and the mean 
separation of absorbers, (d* ep ) = (d sep z1), 24 = (1 + z)/4 
are plotted in Fig. 1 for 1.6 < 2 < 4. These variations are 
well fitted by 

<lgA^> = (\g{N HI /zl)) = 13.3 + 0.08, 

(b) = (26±2.1)km/s, z 4 = (l + z)/4, (10) 

(d* v ) = (2b/H(z)z 3 4 /2 ) = (0.12 ± 0.01)/i _1 Mpc , 

{d* ep ) = {d sep zt} = (1.3±0.15)ft _1 Mpc, 

respectively. For the sample of 7 411 separations we get 

(d* sep ) = (1 ±0.1)fe _1 Mpc. 

Detailed discussion of the observed characteristics of 
absorbers can be found, for example, in Kim, Cristiani & 
D'Odorico (2002), Kim et al. (2002, 2004). 

The observed probability distribution functions, PDFs, 
for the Doppler parameter, P(b) , the hydrogen column den- 
sity, P{Nhi / z\), and the absorbers separation, P{d a e V zl), 
are plotted in Fig. 2. Such choice of variables allows us to 
suppress the redshift evolution of the mean characteristics 
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Figure 1. Rcdshift variations of fraction of measured absorbers, 
fline = ^^abs( z )/^abs> (top panel), mean Doppler parameter, 
(b(z)) and HI column density, (\g(Nnjz^ 2 )) (two middle panels), 
and the mean absorber separation, (d aep z 2 ) (bottom panel). All 
functions are normalized over the mean values given by (10). 

of absorbers. PDFs for so corrected parameters only weakly 
vary with redshift and, as was discussed in Paper I, for the 
main fraction of absorbers, these variations do not exceed 10 
- 15% . These PDFs are well fitted by exponential functions 

Pfit(xHi) ~ 0.4 exp(-0.8x H ) + 2.9 exp(-5x H ) , 

, . J 0.15exp(2.8x 6 ), b < b rap-; /-i -I \ 

mt(Xb) ~ ^ .9cxp(-2.3 a;f) ), b > b rap , 

Pfit(x s ) w 3.5exp(— 1.8a; s )erf 4 (v / l-83; s )/v / x7, 

_ b Nhi/zI _ d aep z\ d* sep 

X "-(&)' XHI ~ {Nhi/zD ' Xs ~ (d sep zl) ~ <d*^> ' 

where again Z4 = (1 + z)/4 and & rap w 23.5km/ s ~ (6t 9 ) 
(3) discriminates between absorbers situated in the increas- 
ing and decreasing parts of the PDF P(xt) in Fig. 2. The 
similarity of & rap and (6& g ) is an independent confirmation 
of estimates (3). For the scatter of measured PDF, P(x a ), 
around the fit (11) we have 

P fit (x s )/P(x s ) « 1.04 ± 0.15, x s < 2 . (12) 

Any model of the forest has problems explaining the 
complex shape of the observed PDF P(xb) and the existence 
of absorbers with b < bb g . These absorbers are sometimes 
related to the unidentified metal lines what is perhaps pos- 
sible for colder absorbers. However, large fraction of such 
absorbers (~ 40%) indicates that majority of them must be 
related to usual hydrogen clouds formed within regions with 
lower background temperature. Our analysis shows that sta- 
tistical characteristics of these absorbers are consistent with 
expected ones for absorbers formed due to adiabatic com- 
pression in rapidly expanded regions (see Sees. 3.5 & 4.5). 
Examples of such absorbers were also found in simulations 
(see, e.g., Bi & Davidsen 1997; Zhang et al. 1998; Dave et 



i 




12 3 

X b> X HI> X s 

Figure 2. Observed PDFs for the Doppler parameter, P(xi,), 
the hydrogen column density, P(xhi), and absorbers separations, 
P(x s ), are plotted together with fits (11) (solid lines). 

al. 1999). This subpopulation can also contain some number 
of "artificial" caustics (McGill 1990) . 

The mean size of rapidly expanded regions, (D rap ), can 
be measured as a distance along the line of sight between 
two absorbers with b > b rap closest to a colder absorber. It 
increases with time as 

(D rap zt) w 1.6 ± 0.2/i _1 Mpc , (13) 
and the typical mass associated with these regions is 

M rap ~ ^(p(z)S m D 3 rap ) ~ lO 13 ^ 3 (5 m )M ^) , (14) 

where 8 m = p m /{Pm) < 1 . At 24 ~ 1 the mass M rap is in 
the range of galactic masses and it rapidly increases with 
time. This result is consistent with the expected symmetry 
of positive and negative initial density perturbations what 
leads to formation of both galaxies and rapidly expanded 
regions. 

Even these results allow one to obtain some inferences 
on the evolution of the forest: 

(i) Regular redshift variations of the mean observed char- 
acteristics of absorbers (10) and weak redshift dependence 
of PDFs (11) indicate the self similar character of absorber's 
evolution over the whole range of redshifts under considera- 
tion. So, we can conclude, that at these redshifts the evolu- 
tion is dominated by a balanced action of the same physical 
factors. 

(ii) The complex form of the PDF P(b) confirms that 
absorbers have been formed within both rapidly and mod- 
erately expanded regions. Analysis of these subpopulations 
of absorbers allows one to estimate some parameters of such 
regions. We will discuss this problem in Sees. 3.5 &4.5 . 

(iii) The very wide range of measured Doppler param- 
eters, 0.2 < b/b bg < 5, indicates a wide variety of initial 
perturbations. Such a wide range of Doppler parameters is 
not reproduced in simulations (Meiksin et al. 2001). 
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(iv) The decline of the measured (Nhi) oc (1 + z) could 
be related to the retained expansion of majority of observed 
absorbers in the transverse directions. 

(v) The growth of the mean observed separation between 
absorbers, (d se p) oc (1 + z)~ 2 indicates the progressive de- 
cline of the number density of observed absorbers that can 
be related to the decrease of their hydrogen column density 
under the observational limit Nhi ~ 10 12 cm~ 2 . 

(vi) Comparison 

of measured {d v (z)} and (1 + z)^ 1 (dse P {z)) (10) indicates 
that at z > 4 the expected overlapping of absorbers be- 
comes essential and, perhaps, absence of observed absorbers 
in spectra of the farthest QSO (Fan et al. 2002, 2003, 2004) 
can be partly related to this effect. 



3 MODEL OF ABSORBERS FORMATION 
AND EVOLUTION 

3.1 Physical model of absorbers 

Many models were proposed during the last twenty years 
to explain regular redshift variations of the mean observed 
characteristics of absorbers (10) and weak redshift depen- 
dence of PDFs (11) (see references in Rauch (1998), and 
in Paper I and Paper II). The simplest one connects the 
absorber characteristics with a suitable set of early formed 
equilibrium clouds. This model naturally explains the weak 
redshift dependence of the mean Doppler parameter, (6), 
and the observed PDFs (11). It explains also the regular 
redshift variations of the mean absorber separation (10). In- 
deed, the mean proper free path between absorbers is 

(1 + zy 1 {d sep (z)) oc (1 + z)~ 3 oc {n ab3 (z)Sabs(z))' 1 , (15) 

where n abs & S abs are the number density and the surface 
area of such clouds orthogonal to the line of sight. In this 
case we have from (15) 

n abs (z) tx (1 + z) s , S ab3 (z) oc const(z) . 

However, it is well known that the existence of such low 
density equilibrium clouds cannot be reasonably explained. 
Moreover, with this model the observed evolution of the 
mean column density of neutral hydrogen, {Nhi} oc (1 + z) 2 , 
must be related to the monotonic evolution of the UV back- 
ground, (Fi2) oc (l + z)~ 2 , what disagrees with the available 
estimates of Ti2. 

Numerous simulations indicate that the majority of the 
LSS elements can be related to the extended anisotropic 
moderate density clouds with a complex internal structure 
and low density envelope. Such long lived clouds are re- 
laxed along the shorter axis and are expanded/compressed 
along transverse directions. The formation of DM pancakes 
as an inevitable first step of evolution of small perturbations 
was firmly established both by theoretical considerations 
(Zel'dovich 1970; Shandarin & Zel'dovich 1989) and numer- 
ical simulations (Shandarin et al. 1995). The anisotropic 
galaxy walls such as the Great Wall are observed in the 
SDSS and 2dF galaxy surveys and represent examples of 
the Zel'dovich pancakes. This suggests that absorbers can 
also be linked with the Zel'dovich pancakes and the reason- 
able model of the forest evolution can be constructed on the 
basis of the Zel'dovich theory. 



However, in some papers devoted to the hydrodynam- 
ical simulations of the forest evolution absorbers are often 
identified with unrelaxed moderate density clouds and their 
Doppler parameter b is related to the gradient of veloc- 
ity of infalling matter rather than to thermal velocities. Of 
course, observed samples include some fraction of such ob- 
jects mainly among weaker absorbers. However, properties 
of richer absorbers with metal systems are found to be in 
a good agreement with a model assuming local hydrostatic 
equilibrium (see, e.g. Carswell, Schaye & Kim 2002; Telferet 
al. 2002; Simcoe, Sargent, Raugh 2002, 2004; Boksenberg, 
Sargent, Raugh 2003; Manning 2002, 2003 a,b; Bergeron 
& Herbert-Fort 2005). This means that, possibly, the frac- 
tion of simulated extended unrelaxed clouds is artificially 
enhanced by technical limitations (see, e.g., discussion in 
Meiksin, Bryan & Machacek 2001; Manning 2003 and Paper 
II). It is important that these simulations cannot yet repro- 
duce well enough the characteristics of observed absorbers 
discussed in previous Section. 

This discussion shows that an adequate physical model 
of the complex evolution of the forest has not been proposed 
yet. Indeed, the relaxation of DM pancakes leads to com- 
plex internal structure of pancakes, the adiabatic compres- 
sion and/or expansion of absorbers in transverse directions 
is changing their overdensity and temperature, the radia- 
tive cooling and bulk heating leads to the drift of the gas 
entropy and overdensity but leaves unchanged the depth of 
the potential well formed by DM distribution. Merging of 
pancakes increases more strongly the depth of the potential 
well and the gas entropy but the overdensity of the gaseous 
component increases only moderately. The temperature and 
overdensity of the trapped gas are rearranged in accordance 
with the condition of hydrodynamic equilibrium across the 
pancake continuously. These processes imply the existence 
of a complicated time- dependent internal structure of ab- 
sorbers. 

In this paper we discuss the model of absorbers evolu- 
tion based on the Zel'dvich theory. This analytical model 
links the self similar evolution of dominated DM component 
with observed characteristics of absorbers. Here we assume 
that: 

(i) The DM distribution forms a set of sheet-like clouds 
(Zel'dovich pancakes), their basic parameters are approxi- 
mately described by the Zel'dovich theory of gravitational 
instability applied to the CDM or WDM initial power spec- 
trum (DD99; DD04). The majority of DM pancakes are 
partly relaxed, long-lived, and their properties vary due to 
the successive accretion of matter, merging and expansion 
and/or compression in the transverse directions. 

(ii) Gas is trapped in the gravitational potential wells 
formed by the DM distribution. For majority of absorbers, 
the observed Doppler parameter, 6, traces the gas tempera- 
ture and the depth of the DM potential wells. We consider 
the possible macroscopic motions within pancakes as sub- 
sonic and assume that they cannot essentially distort the 
measured Doppler parameter. 

(iii) The gas is ionized by the UV background and for the 
majority of absorbers ionization equilibrium is assumed. 

(iv) For a given temperature, the gas density within the 
potential wells is determined by the gas entropy created 
during the previous evolution. The gas entropy is changing, 
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mainly, due to relaxation of the compressed matter, shock 
heating in the course of merging of pancakes, bulk heating 
by the UV background and local sources and due to radia- 
tive cooling. These processes slowly change the entropy and 
density of the trapped gas. Random variations of the inten- 
sity and spectrum of the UV background enhance random 
scatter of the observed properties of absorbers. 

In this model we expect that the mean surface density 
of both DM and baryonic components of pancakes weakly 
varies with time and the mean density of matter compressed 
within absorbers decreases as (p a bs) oc (1 + z) v , v ~ 1.5 — 1 
for colder and hotter absorbers, respectively. We cannot dis- 
criminate between evolution of the absorbers number den- 
sity, n a bs, and their surface area, S a bs- However, numerical 
simulations demonstrate that the general tendency of evo- 
lution is a sequential growth of masses and sizes of pan- 
cakes accompanied by a drift of weaker ones under the ob- 
servational limit and decrease of the comoving number den- 
sity of observed absorbers. Thus, as is seen from (15), if 
the surface area of absorbers increases {S a b s ) oc (1 + z)~ K 
then the number density of observed absorbers decreases as 
(n a bs) oc (1 + z) 3+K . We show that the observed evolution 
of mean hydrogen column density and mean absorber sepa- 
rations (10 & 11) coincides with theoretical expectations for 
Gaussian initial perturbations. 

The wide range of the observed Doppler parameters, 
5-10 km/s < b < 100 km/s, demonstrates a complex com- 
position of the forest. In this paper we roughly identify three 
subpopulations of absorbers, namely, hot absorbers formed 
in the course of merging and shock compression, warm ab- 
sorbers formed due to adiabatic and weak shock compres- 
sion in moderately expanded regions, and an unexpectedly 
rich subpopulation of colder absorbers. We link the mea- 
sured Doppler parameter with the depth of ID potential well 
formed by the compressed DM component and we neglect 
the contribution of macroscopic velocities what restricts the 
achieved precision of our approach. Non the less this model 
reproduces quite well the self similar evolution of absorbers 
and allows one to reconstruct the initial correlation func- 
tion of velocities and the initial power spectrum down to 
very small scales. 



3.2 The initial power spectrum and correlation 
functions of the initial velocity field 

In this Section we summarize the main results obtained in 
DD99 and DD04 concerning the evolution of DM pancakes 
and in Sec. 3.9 we show how to improve the estimates of 
the correlation function of the initial velocity field, and the 
shape of initial power spectrum discussed in Paper II. 

As a reference power spectrum of initial perturbations 
we take the standard CDM-like spectrum with the Harrison 
- Zel'dovich asymptotic, 



A 2 k , 



P ( k ) = —T 2 (v)Dw(v), 



V 



k u 
ko 



Mpc 



(16) 



where A is the dimensionless amplitude of perturbations, k is 
the comoving wave number. The transfer function, T(rj), and 
the damping factor, Dw{t)), describing the free streaming 
of DM particles were given in Bardeen et al. (1986). For 



WDM particles the dimensionless damping scale, Rf, and 
the damping factor, Dw , are 

i f n m h 2 kev \ 4/3 2 

R f = 5 I Mdm I , D w = ex.p[-r)R f - (r)R f ) } , 

where Mdm is the mass of WDM particles in keV (Bardeen 
et al. 1986). This relation illustrates the clear dependence of 
characteristics of small scale perturbations on the mass of 
DM particles. 

For the spectrum (16) the coherent lengths of velocity 
and density fields, l v & l p , are expressed through the spectral 
moments, m_2&mo, (DD99) as follows: 



lv — 



1 



ko Q. m h 2 



f ' ° Mpc « 31.4/1" 1 Mpc, (17) 



m_2 



m 



dr, V T 2 (r))D w ( V ) « 0.023, l p = q L , 



p oo 

/ dr, rj 
Jo 



3 T 2 ( V )D W (V), qo = $ — 
mo 



The moment mo depends upon the mass of dominant com- 
ponent of DM particles and as is shown in Sec. 4.2, go < 
10~ 3 , m > 2.4, M DM > 1 MeV . 

As was demonstrated in DD99 and DD04, the basic 
statistical characteristics of structure are expressed through 
the normalized longitudinal correlation function of the ini- 
tial velocity field, v(q) , 

g.(i. g )=3 <(q - vfe » ( r Vfe))> . q=^7^. (18) 



lv 



Here qi & q2 are the unperturbed coordinates of two parti- 
cles at z = 0, q = |q|, and a 2 is the velocity variance. This 
function is expressed through the power spectrum: 

Zv{l v q) = —{ dr) r) 2 cos x [ %T 2 {y)Dw(y), (19) 

m -2 .l j ri y 



V T 2 ( V )D w (ri) 



V m -2 



f 

Jo 



(2 cos a; + x sin x)£ v (l v q)dq , 



x — kl v q, rj = k/k , £„(0) = 1, 



f 

Jo 



dq £{l v q) = . 



Similar relations can be also written for any initial power 
spectrum P(k) . 

For the CDM - like spectrum (16) with q Q < 10~ 3 and 
for the most interesting range 0.5 > q the velocity correla- 
tion function can be fitted as follows: 



£v{q) = (,cdm ~ 1 



1.5q 2 



^2.25q^+q 2 +pl A q°- e +q 2 o ' 



(20) 



where p w 1.1 • 10 — 2 and q was introduced in (17). Further 
on, we will use this function as the reference one and will 
compare it with observational estimates of £«(<?). 

As was shown in DD99 and DD04, the main DM char- 
acteristics depend upon the self similar variable 

2 



q 



4r»(z)[l- £„(«)]' 



(21) 



where the 'time' t(z) describes the growth of perturbations 
due to the gravitational instability. For the ACDM cosmo- 
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logical model (1), and for z > 2 we have 
. / l-MAA 1 / 3 1 1.27r 

r(2) " r(, l^^rj TTl TTT ' (22) 

and to characterizes the amplitude of initial perturbations. 
It is proportional to as, the variance of the mass within a 
randomly placed sphere of radius 8h~ Mpc. Latest estimates 
(Spergel et al. 2003, 2006; Viel et al. 2004b) for the model 
(1) are 

&8 



(j 8 « 0.9 ± 0.1, T = 0.21o-8 w (0.19 ± 0.02) 



0.9 



(23) 



Here 5 m is the mean overdensity of compressed matter above 
the background density. 

For the transverse size of absorbers, the expected char- 
acteristics were estimated in DD04 as follows: 



qtr 



(27) 



Ntr(Ctr) « -^= exp(-Ci 2 r), Ctr ■ - 

(Ct) « 1/2, (J„gtr) « 3r 2 (2)Z„ « 0.45zJ 2 ft _1 Mpc . 

Here again 24 = (1 + z)/4, 2 is the observed redshift and 
l v q tr is the expected size of absorbers at 2 = 0. 



3.3 Expected characteristics of DM absorbers 

In this section we introduce (without proofs) the basic char- 
acteristics of DM pancakes as a basis for further analysis. 
For more details see DD99 and DD04. 



3. 3. 1 DM column density of absorbers 

The fundamental characteristic of DM pancakes is the di- 
mensional, fj,, or the dimensionless, q, Lagrangian thickness 
(the dimensionless DM column density) : 



(p m {z))l v q _ 3H { 



(1 + *) 



p=l v il m (l + zfq , 

7TO 



(24) 



where l v is the coherent length of initial velocity field (17). 
The Lagrangian thickness of a pancake, l v q, is defined as 
the unperturbed distance at redshift 2 = between DM 
particles bounding the pancake (18). 

As was found in DD99 and DD04 for Gaussian initial 
perturbations, the expected probability distribution func- 
tion for the DM column density is 



N g (C) 



W q {< C)=erf 2 (C) ; 



(25) 



«> « \ + ~ « 0.82, (q) a 6r 2 (C) 

Z TV 



2±O2 10 _ 2 



where 24 = (1 + z)/4, (,(q,z) was introduced in (21) and 
W q (< () is the cumulative probability function. 

Strictly speaking, the relations (24) and (25) are valid 
for pancakes formed and observed at the same redshift 
Zots = Zf because after pancake formation the transverse 
expansion and compression changes its DM column density 
and other characteristics. However, owing to the symme- 
try of moderate distortions of general expansion in trans- 
verse directions, these processes do not change the statis- 
tical characteristics for majority of pancakes observed at 
redshift z obs < Zf (DD04). This means that statistically 
we can consider each pancake as created at the observed 
redshift. However, the symmetry is distorted for rapidly ex- 
panded regions observed as absorbers with 6 <C b bg because 
the cross-section of strongly compressed absorbers is small 
and, therefore, they are rarely observed. More details are 
given in Sees. 3.5, 3.6 and in DD04. 

3. 3. 2 Proper sizes of absorbers 

The actual thickness of a DM pancake is estimated as 



(1 + z)8 r 



-, 5m — pm/ {pm(z)) . 



(26) 



3.3.3 Absorbers separation 

An important characteristic of the distribution of absorbers 
is their separation determined as the distance along the line 
of sight between centers of neighboring absorbers observed 
with the separation Az, 



5.5 • 10 J 



Az 




1 Mpc. 



(28) 



(1 + 2) 3 / 2 

This separation is identical to the free path between ab 
sorbers and is quite similar to the popular descriptor 

AN ab3 



Az 



oc {d sep (z))~ 



(29) 



For Gaussian initial perturbations the PDF for the ab- 
sorber separations in Lagrangian space, g!l, is 



N l (Cl) w 4tt- 1/2 exp(-Ci)[(l + 2C L )I>»(Va 

w 2.82exp(-Ci>rf 4 (CL)/v / a , 

Cl = C(Ql,z), q L = d L /l v , (Ci)wl.5, 

where ((q, 2) was introduced by (21) and 



(30) 



D w (x) 



dyexp(y 2 - x 2 ) 



is the Dawson function. 

Comparing the pancake surface density, /n(z), with the 
absorbers separation, d sep , we can also estimate the fraction 
of matter accumulated by absorbers as 



^ p.(l + z) _ l v q ^ (q) 

Jabs — , » , , 7 \Jabs) ~ , , 

\pm)(lsep dsep W^/ 



0.55. 



(31) 



However, this formal estimate is of limited significance be- 
cause it does not consider the complex processes of absorbers 
evolution. More detailed 3D analysis (DD04) shows that at 
a later period of absorbers evolution we can expect 



f abs ~ 0.38 - 0.44 



(32) 



The observed characteristics of absorbers are measured 
in the redshift space where both the proper motions of ab- 
sorbers and their peculiar velocities distort the PDF (30) 
and, in particular, lead to the merging and artificial blend- 
ing of absorbers. These processes are driven by the spatial 
modulations of gravitational potential formed by the large 
scale perturbations and can be also described in the frame- 
work of Zel'dovich approximation (DD99). However, such 
description depends upon the size of absorbers (26) which 
variations with redshift cannot be described theoretically. 
So, characteristics of absorbers in both real and redshift 
spaces can be determined only approximately. 
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For the PDF of absorbers separation in the redshift 
space we get 

N sep « 2.4cxp(-1.35x rd )[l - 0.85 exp(-1.35a; rc j)] (33) 

_ CiQrd, Z) iff \\ nil dsep 

x rd = jf {(,{qrd,z)) « 2.13, q rd =—^, 

(d sep ) « 6/„r 2 (C(grd,«)) « 1.3«r 2 ft _1 Mpc(T /0.2) 2 , 

where z 4 = (1 + z)/4 and £(<?, z) & d sep are defined by (21) 
and (28). These expressions describe correctly properties of 
larger separations but become unreliable for smaller separa- 
tions where the influence of the proper sizes of absorbers is 
more important. In the real space the expression (33) also 
approximates the PDF of distances between neighboring ab- 
sorbers but the mean value {(real) ~ 1-64 is smaller than 
that in the redshift space. 

The same approach allows us to determine the expected 
characteristics of merged absorbers. Instead of (25), in red- 
shift space, we get for such absorbers 

N mrg (x) « 1.3exp(-l.lx)[l-0.5exp(-2.2a:)], (34) 

X — C,mrg/ (C mr 9*) s C rnr 9 — C(^? mr ff> ^Oi (C mr <?) — 1.22 , 

where ((q,z) was introduced in (21) and q mrg is the di- 
mensionless DM column density of merged absorbers. This 
relation indicates that the population of poorer merged ab- 
sorbers is suppressed and, in particular, the mean DM col- 
umn density of merged absorbers is larger than that for all 
absorbers as given by (25). 

The relations (25), (31), (33), (34) show that during the 
self similar period of structure evolution, when the relations 
(20) and (22) are valid, we can expect regular variations of 
the basic characteristics of absorbers such as their DM col- 
umn density, q, separation, d aep , and fraction of matter accu- 
mulated by absorbers, f abs . These regular variations are dis- 
torted at small redshifts when the growth of perturbations 
is decelerated and at higher redshifts, when the blending of 
absorbers becomes more important. 



3.4 Doppler parameters of absorbers 

For relaxed and gravitationally confined absorbers their 
Doppler parameters are closely linked to the potential wells 
formed by the DM distribution. As is well known, for an 
equilibrium slab of DM the depth of its potential well is 



A<E> : 



wo 



(p(z))S, 



-04 



3 2 £ 2 (l + z) n 
-v -z 6$ . 



(35) 



HqI v \ 



1 720km/ s 



0.15 

h 2 



where the random factor 0<i> characterizes the inhomogene- 
ity of DM distribution across the slab and the evaporation 
of matter in the course of its relaxation. 

Analysis of numerical simulations (Demiahski et al. 
2000) indicates that the relaxed distribution of DM com- 
ponent can be approximately described by the polytropic 
equation of state with the power index 7 m « 1.5 - 2 . Thus, 
for 7 m = 2 the equilibrium density profile across the slab 
can be directly found and 

9* = 4/tt w 1.3 . 



The actual distribution of DM component across a slab and 
the value of O* depends upon the relaxation process which is 
essentially accelerated by the process of pancake disruption 
into the system of high density clouds and filaments. In the 
course of relaxation ~ 10 - 15% of matter is evaporated 
what decreases the factor 0$. This means that the expected 
6$ ~ 1 randomly varies from absorber to absorber. 

The Doppler parameter is defined by the depth of po- 
tential well (35) and for the isentropic gas with p gas oc p% 3 s 
trapped within the well, we get 



bl 3 + ^ 2 *{l + z)Q, 



5~ -» ' 10-°* m v - ' (36) 

Variations of the gas entropy across absorbers increase the 
random variations of & b. 

For hot absorbers with 6 ^> (b bg ) we can neglect the 
difference between the actual and mean background tem- 
perature and in (36) use (b bg ) instead of b bg . In such a way 
we link q 2 /8 m with bk, {b bg ) with a reasonable precision. 



3.5 Absorbers within rapidly expanded regions 

For significant fraction of absorbers - up to 20% - the 
Doppler parameter is smaller then the expected mean back- 
ground one (3). Such absorbers are often related to unrecog- 
nized metal lines. However, both theoretical arguments and 
numerical simulations show that such absorbers can also be 
related to hydrogen clouds formed within colder rapidly ex- 
panded regions. 

The temperature of relaxed HI absorbers formed by 
the compression of matter (36) cannot be smaller than the 
background one (3). However, within low density rapidly ex- 
panded regions the background temperature given by (5) is 
smaller than the mean one, and in such regions the hydro- 
gen clouds with b bg < b < (b bg ) can be formed. Our analysis 
(Sec. 4.5) shows that majority of colder absorbers could be 
related to such hydrogen clouds. 

According to the Zel'dovich theory of gravitational in- 
stability and for Gaussian initial perturbations, regions with 
moderate distortions of the cosmological expansion domi- 
nate and probability to find rapidly expanded or compressed 
regions is exponentially small. However, the observations of 
galaxies at redshifts under consideration corroborate the ex- 
istence of strong distortions of cosmological expansion at 
least on galactic scales. This means that owing to the sym- 
metry of positive and negative initial density perturbations 
we should also observe rapidly expanded regions with small 
background density and temperature. Number of such re- 
gions exponentially decreases for larger distortions of the 
expansion. 

Owing to the same symmetry, for majority of absorbers, 
fraction of absorbers with a moderate random expansion 
and compression in the transverse directions are close to 
each other and the influence of this factor weakly distorts 
the mean absorbers characteristics. However, within rapidly 
expanded regions all absorbers are adiabatically expanded 
in the transverse directions what, in particular, decreases 
the fraction of neutral hydrogen, 

x H oc n b /b'' /2 T 12 oc nl /2 /r 12 , 

and its observed column density, and leads to a systematic 
drift of absorbers under the observational limit IgNui ~ 12. 
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Theoretical estimates show that owing to the correla- 
tion of velocity perturbations across the pancake and in the 
transverse directions the rate of pancakes formation within 
rapidly expanded regions is smaller than the mean one. For 
the CDM like initial power spectrum the coefficient of cor- 
relation of orthogonal velocities is c v « 1/3 (DD04). In this 
case, for the fraction of matter, f rap , and the mean column 
density of DM component, q rap in the rapidly expanded re- 
gions, we expect 



0.33(/ Oi .>, ((rap) = (((qra P , z)) » 0.5(C) ■ 



(37) 



3.6 Characteristics of the gaseous component 

The observed column density of the neutral hydrogen can 
be written as an integral over the line of sight through a 
pancake 



Nhi 



-I 



dx p b x H = 2(xh) 



(n b (z))l v q 0.5 

1+2 COSO 



(38) 



Here (xh) is the mean fraction of the neutral hydrogen and 
cosO takes into account the random orientation of absorbers 
and the line of sight ((cos8) « 0.5). As was noted in Sec. 
3.4, we assume also that both DM and gaseous components 
are compressed together and, so, the column densities of 
baryons and DM component are proportional to each other. 

Under the assumption of ionization equilibrium of the 
gas (6) and neglecting a possible contribution of macroscopic 
motions to the fc-parameter (T oc b 2 ), for the fraction of 
neutral hydrogen and its column density we get: 



(x„) = x 8 b p-' i,2 {l + «r i4 6„ 13 = b/(b bg ) , 



33/14, 



qSb 



No 



ri 2 /3 3 /- 



■,61/14 



No 



5.5 ■ 10 12 cm 2 Qh, 



(39) 



the thickness of a gaseous pancake is estimated as 
l v q 



da 



(l + z)S b 



(41) 



For adiabatically compressed absorbers it is larger then the 
thickness of the DM pancake but for shock compressed ab- 
sorbers they are close to each other. So defined d abs can be 
compared with the estimated redshift thickness of absorbers 
determined by the observed Doppler parameter, 



d v = 2b/H(z) . 



(42) 



We can expect that, as usual, d v > d abs owing to the impact 
of thermal velocities. In spite of the limited precision of de- 
termination of d v & dabs these estimates allow us to correct 
the derived DM column density, q. 

For long lived absorbers the influence of the bulk heat- 
ing can be estimated in the same manner as it was done for 
characteristics of the background (Demianski & Doroshke- 
vich 2004b). Solving the equation of thermal balance for 
absorbers formed at z = Zf and observed at z « z; wc 
obtain for the entropy of compressed gas: 



F 3/2 
bg 



7/3V2(2) 



1 + z 

l + Zf 



3/2 



(43) 



As is seen from this relation, the bulk heating is negligible for 
shock compressed absorbers with F s (zf) 3> F bg (z), /3(z) 3> 
1. For adiabatically compressed long lived absorbers with 
z f > 2, F s (z f ) « F bg (z f ) < F bg (z) we get 



.3/2, 



bg 



1/2, 



7/2 



(44) 



This result demonstrates that the bulk heating is specially 
important for absorbers with j3< 1 formed within rapidly 
expanded regions. For such absorbers the merging is sup- 
pressed and their number decreases mainly owing to the 
drift under the observational limit. 



0i 



9s 0.15 (cose) ( n b h 



Q 3 J 7 Qr, 
bg 



h 2 cosd 



0.02 



8 b = n b /(n b ) . 



where Fi2, (b bg ), Q bg and xo were defined in (3) and (7) and 
the factor O x ~ 1 describes the inhomogeneous distribution 
of ionized hydrogen along the line of sight. 

However, the overdensity of the baryonic component, 
6 b is not identical to the overdensity of DM component, 6 m , 
(see, e.g., discussion in Matarrese & Mohayaee 2002). In- 
deed, the gas temperature and the Doppler parameter are 
mainly determined by the characteristics of DM component 
(36) but the gas overdensity is smaller than that of DM com- 
ponent due to larger entropy of the gas. Moreover, the bulk 
heating and cooling change the density and entropy of the 
gas trapped within the DM potential well. These processes 
change the baryonic density of pancakes and we can write 



6 b = Q b (z)6 m , e b (z)<l. 



(40) 



The factor <d b should be small for absorbers formed due to 
adiabatic and weak shock compression because of the large 
difference between entropies of the background DM and the 
gas, and O b — > 1 for richer hot absorbers formed due to 
strong shock compression when entropies of both compo- 
nents are comparable. 

Similarly to the proper thickness of a DM pancake (26), 



3.7 Observed characteristics of absorbers 

Eqs. (36) and (39) relate three independent variables, 
namely, q,6 m &i6 b . To find the DM column density, q, it 
is therefore necessary to use an additional relation which 
connects the basic parameters of absorbers. Here we assume 
that the richer absorbers with b > b t hr = I3thr(b bg ) > (b bg ) 
are formed due to shock compression and this process is ac- 
companied by strong relaxation of compressed matter. For 
such absorbers 6 b ~ 6 m and their DM column density is: 



6 b : 



Vo(l + z) 2 



b 2 " (bbg) 2 



q , b > bthr 



N H ir 12 b 2 - (b bg ) 2 ( b 



N 



\(bbg) 



3/2 



(l + z)- 75/14 



(45) 



Formation of absorbers with b < b t hr = Pthr(b bg ) 
is accompanied by adiabatic or weak shock compression 
of baryonic component. Assuming that the compression of 
baryons is described by the polytropic equation of state with 
76 = 5/3, we can expect that for recently formed absorbers 
with b< b thr , F s (z) w (F bg (z)), 

N H iT 12 



/3- 3/2 (l + 2) 



No 



-61/14 



(46) 
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and for long lived absorbers (44) 



iVo 



p- 2 (l + z) 



-61/14 



7/2 



(47) 



The relations (45) - (47) determine the dimensionless 
column density of DM component corrected for the impact 
of gaseous pressure. These relations can be successfully ap- 
plied to absorbers formed due to adiabatic and strong shock 
compression with various degrees of relaxation. However, the 
boundary between these limiting cases must be established 
a priory. So, to discriminate absorbers described by (45) 
and (46, 47), we use the threshold Doppler parameter, bthr 
which, in fact, characterizes the Mach number of the inflow- 
ing matter. Thus, absorbers with b < b t hr, 



bthr — Pthr(b bg (z)), (3thr « 1-5 — 2 , 



(48) 



can be conveniently considered as the adiabatically com- 
pressed while absorbers with b > bthr arc considered as 
formed by shock compression and they contain strongly re- 
laxed matter. 

The precision of these estimates is moderate and the 
some uncertainties are generated by the poorly known T 7 
and the parameters 0$ and Oh, which vary - randomly 
and systematically - from absorber to absorber. Estimates 
of these uncertainties can be obtained from the analysis of 
the derived absorbers' characteristics. However, the main 
uncertainties in the estimates of q come from the unknown 
cos 9, and, for strongly relaxed shock compressed absorbers, 
from modulations of bt g . 

Indeed, for adiabatically compressed absorbers q oc cosO 
and Eqs.(46, 47) underestimate q for cos# > (cos 9) = 0.5 
and overestimate it for cos 9 < (cos 9) . To reveal and to cor- 
rect the most serious uncertainties we will use the condition 
dabs/d v < 1 where the real, d a bs, and the redshift, d v , size 
of absorbers were introduced in (41, 42). For absorbers with 
dabs > d v we will substitute the 'true' column density, q t 
denned by relation 



q t = q d v /d a bs, d a 



d v , 



(49) 



instead of the measured one (46) or (47) . Perhaps, more de- 
tailed reconstruction of the shape of absorption lines could 
allow us to reveal stronger deviations from the Doppler pro- 
file and, so, to identify the influence of absorbers orientation 
and macroscopic velocities. 

For shock compressed absorbers the influence of cos 9 is 
not so strong as q oc cos 1 / 3 6 and the expected modulation 
of bb g is more important. Indeed, these absorbers are formed 
owing to the merging which is more probable in slowly ex- 
panded regions with bb g > (bbg)- Thus, for absorbers with 
b > bthr, d abs > d v , we will also correct q&id abs by relation 
(49). 

These corrections essentially improve estimates of the 
correlation function of initial velocity field. 



3.8 Regular and random variations of absorbers' 
characteristics 

The most fundamental characteristic of absorbers is their 
DM column density, £(q, z) oc q(l + z) 2 (21). It depends on 
the process of formation and merging of pancakes, is only 
weakly sensitive to the action of random factors and defines 
the regular redshift variations of absorbers characteristics. 



For shock compressed absorbers, the evolutionary his- 
tory of each pancake and the action of random factors dis- 
cussed in the previous subsections are integrated in the en- 
tropy of the baryonic component, 



S b = lnF.(z) = S bg + 2/3ln(l3 3 /5 b ) . 



(50) 



If the structure of a relaxed DM pancake can be described 
by the polytropic equation of state with the effective power 
index 7 m then we can introduce also the entropy of DM 
component, S m - For probable value ) m « 1.5 - 2, entropies 
of DM and baryonic components are quite similar to each 
other, S m ~ Sb- This means that for the strongly relaxed 
shock compressed absorbers the evolutionary history is char- 
acterized quite well by two functions, £ and Sb- 

For adiabatically compressed absorbers the baryonic en- 
tropy is identical to the background one given by (4) while 
the observed 6 and Nhi depend upon the distribution of the 
compressed DM component. For such absorbers the PDFs 
and the random scatter of observed characteristics are de- 
fined mainly by variations of the expansion rate and the 
background density and temperature. These characteristics 
as well as the entropy and overdensity of compressed DM 
component, S m &<5 m , now cannot be derived from observa- 
tional data with a reasonable reliability. This problem de- 
serves further investigations. 



3.9 Reconstruction of the initial power spectrum 

The basic relation of Zel'dovich theory of gravitational in- 
stability can be suitably written for the difference of coor- 
dinates of two particles, 



An 



lyT(z) 

1 + z 



[ qi /T( Z )-ASi(q)/lv], » = 1,2,3, 



(51) 



where n are the Euler coordinates of particles, l v , t(z) and 
qi were introduced in (17), (18) & (22) and Si(q) is a random 
displacement of a particle with respect to its unperturbed 
position. As is seen from this relation, for pancakes AS w 
l v q/r(z) and, so, some statistical characteristics of AS can 
be obtained by measuring q/r(z). 

Using the measured redshift, z, and DM column den- 
sity of absorbers, q, we determine the cumulative PDF of 
absorbers W (, s [> q/r(z)] and, for each q/r(z), we compute 
(q) and a 2 = (q 2 ) — (q) 2 . For a chosen W q (C) (25), we solve 
numerically the equation 

W ob s[< q/r(z)] = W q (0 = erf 2 (C) , (52) 
with respect to C(l, T ) and, thus, we obtain the function 



(53) 



4r 2 C(g) ' 

For the most interesting range q/r(z) <C 1,(< 1 we have 

W 9 (C)«4C77r, 1 - oc q 2 /y/W q (C) ■ (54) 

The same approach can be applied to the absorbers 
separation, d sep , and £ r£ j = C(<jw, z), &iq r d introduced by 
(33) . In this case instead of (52) and (54) we have 

W ob s[< d sep /lvT(z)] = W S ep(< (rd) ~ 

1 - 1.778 exp(-0.634Crd)[l - 0.425 exp(-0.634Crd)] , (55) 
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2 
Ird 



1 " iv{qrd) '- 4^M' 

and, for q rd /T(z) < l,( rd < 1, we get 

W/ sep (Crd) « 0.17Crd, 1 - £v OC Ijrd/Wsep ■ 

The observed functions P bs(q) and Pobs(d 3e p) can be 
compared with corresponding expectations (25) and (33). 
The correlation function £„(<?) derived from (53) and (55) 
can be compared with the reference function £,cDAi(q) (20). 



4 MODEL DEPENDENT STATISTICAL 
CHARACTERISTICS OF ABSORBERS 

In this Section, evolution of the basic model dependent char- 
acteristics of absorbers is discussed in the ACDM cosmologi- 
cal model (1), and for the background temperature, T bg , the 
Doppler parameter, b bg , and the entropy function F bg given 
by (3) and (4) . 

4.1 Parameters of the model 

The model of absorbers discussed in Sec. 3 includes poorly 
known random parameters, Bbg, &<s>&zB x , which cannot be 
estimated a priory, this leads to a moderate random scat- 
ter of derived absorbers characteristics. Further on we will 
assume that 



B b0 = 6$ = e a 



i . 



(56) 



However, as was discussed in Sec. 3.7, the main sources of 
uncertainty in (45 - 47) are the random orientation of ab- 
sorbers with respect to the line of sight, measured by cos 6 
and, for strongly relaxed shock compressed absorbers, ran- 
dom spatial modulation of the expansion rate and the back- 
ground temperature, b bg > (b bg ). Large distortions of the 
derived parameters of absorbers can be revealed and partly 
corrected using the relation (49). However moderate distor- 
tions of the same parameters that can not be easily corrected 
restrict the precision of our approach. 

The poorly known radiative ionization rate, Ti2, (8) is 
also an important source of uncertainty. As was noted in Sec. 
2.1, now there are approximate estimates of the UV back- 
ground produced by the observed QSOs but the expected 
ionization rate should be corrected for the absorption and 
reemission of UV radiation by the gas compressed within 
high density clouds (Haardt & Madau 1996) and for the 
additional emission of UV radiations by galaxies (at z < 2) 
and at z > 3 by poorly known sources such as Ly-a emitters 
(Boksenberg, Sargent & Rauch 2003; Giavalisco et al. 2004; 
Ouchi 2005). 

Below we will describe the ionization rate by the ex- 
pression 



Go = 4.3, 



exp 



2o-2 



(57) 



1.0, 



= 1.58, p 7 = 1.5, 



where the choice of p 7 ,z 7 and er 7 corrects T 12 (z) for the 
impact of additional sources of radiation. 

Relatively small value of Ti2(2) « 3 (57) is close 
to the observational estimates of Scott et al. (2002) and 
Ti2(5) « 0.2 is similar to the estimates of Fan et al. (2002, 



2004). It shows that at z ~ 2 the UV background was proba- 
bly overestimated in Haardt & Madau (1996) and Demianski 
& Doroshkevich (2004b) . Indeed, relatively small observed 
(Nhi) ~ 10 13,4 cm~ 2 , with lg Nhi < 15, and moderate num- 
ber of lines Nu ne — 226 with IgNni > 15 in our sam- 
ple, shows that the absorption of the UV background by 
HI cannot be very important. However, stronger absorp- 
tion by Hell found by Levshakov et al. (2003) can dis- 
tort the spectrum of the UV background and decrease Go 
and the background temperature, T bg & b bg (3). Below the fit 
(57) will be tested by comparing the theoretically expected 
and derived functions, £„(g) and £ v (q s ), (C( 2 )) ~ 0.82 and 
(l v q(z)) «0.44(d sep («)}. 

For such Ti2(z) we get for the Gunn-Peterson optical 
depth, T G p(z), 



top (2) a 0.045/L™, t gp (5) » 5/L 



(58) 



where fhom is the fraction of homogeneously distributed 
matter. Estimates obtained in Sec. 4.3 show that at z < 
4, fhom ~ 0.5. However, at such redshifts tgp is measured 
for regions with suppressed lines where larger Pi2 can be 
expected. For large z > 4 we can take fhom ~ 1 with large 
scatter. 

The choice of Ti2 (z) (57) coincides with 
T ^0.19, o- 8 -0.9, 

what agrees quite well with a$ derived by Spergel et al. 
(2003, 2006) and with the independent estimate (10). The 
final results only weakly depend upon the threshold param- 
eter bthr ~ (1.5 — 2)b bg discriminating between adiabatically 
and shock compressed absorbers. Here we use bthr ~ l-5b bg . 

Such choice of the model parameters allow us to obtain 
reasonable description of properties of absorbers for the se- 
lected sample. Variations of limits used for the sample se- 
lection lead to moderate variations of parameters (57) . 

4.2 Correlation functions of the initial velocity 
field 

4-2.1 Correlation functions derived from absorber 
separations 

Using the method described in Sec. 3.9 and characteristics 
of separations between absorbers obtained in Sec. 2.2 we can 
estimate also the correlation function of the initial velocity 
field, £ v (q r d) ■ This approach uses only the measured red- 
shifts and, so, the derived function does not depend upon 
the measured b & Nhi and the model of absorbers discussed 
in Sec. 2. However, the theoretically derived PDFs (33) and 
(55) become unreliable at small separations and, as was 
noted in Sec. 2.2, already at redshifts z ~ 3 the blending of 
absorbers becomes essential what distorts the derived func- 
tion £ v (q r d)- 

The observed cumulative PDF, in redshift space, 
W bs[qrd/T~(z)], and the reconstructed correlation function 
1 — (,v(q r d) (55) are plotted in Fig. 3 for the samples of 19 
QSOs (6 251 and 7411 separations) and 14 QSOs (3 660 sep- 
arations). For these samples we have, respectively, 

W obs [< qrd/T(z)]/W sep (< Ou) « 1-07 ± 0.19, 1.03 ± 0.2 , 



- £,v{q r d) 
■ tfit(qrd) 



0.95 ±0.25, 1±0.2, 



(59) 
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Figure 3. Top panel: the cumulative PDF, W t a (q r d/ T ), together 
with the fit Wsep(Crd) (55) (dashed line). Bottom panel: the cor- 
relation function of initial velocity field, 1— £v(Qrd)i together with 
the theoretical fits (20) (dashed line) and the fit (59) (solid line). 
Points and stars show the same functions found for the sample of 
14 QSOs and for the sample of 7 411 separations in 19 spectra, 
respectively. 



■ £/«(<7rd) 



y/2.2S<£ d + <£ d +p\-* q ™- 



2.1 • 10" 



0.5 > <? > 1.4 • 10 3 , 17ft *Mpc > l v q r d > 0.03ft : Mpc. 

where C/« = ({£,fit,r) is given by (21). 

As is seen from Fig. 3, at q r d > 10~ 2 , l v q r d > 
0.33ft _1 Mpc the derived correlation function coincides with 
the standard CDM - like one (20) for heavy DM particles, 
with go < 10~ 3 , Mom > lMeV . Similar results can be found 
also with the PDF (30) what indicates the weak sensitivity 
of the derived correlation function Q it on the detailed shape 
of the used PDF. 

At small scales, 10 -3 < q r d < 10~ 2 , we see an excess of 
power with respect to the reference function (20). This 
excess is caused by the deficit of observed absorbers with 
small separation which increases the derived function, 1 — 
£,v(q r d) oc W~Jp (see Sec. 3.9). In this range of q r d the number 
of measured separations is limited, N 3 e P ~ 320 and N 3 e P ~ 
170 for the samples of 19 and 14 QSOs, respectively, what 
decreases reliability of this result. To test impact of this 
factor we calculated the function £ v (q r d) for the sample of 
7411 separations obtained without placing any restrictions 
on the properties of absorbers. For this extended sample 
the number of small separations increases up to 630 but 
the difference between the derived and reference functions, 
£,v(q r d) and £cDM(q r d), at q r d < 10~ 2 remains the same. 

These results indicate that possibilities and applicabil- 
ity of this approach are limited. Indeed, the positions of 
absorbers are measured in the redshift space with a typical 
error Az ~ 5-10" 5 , Aq rd ~ 3- 10~ 2 [4/(l±z)] 3/2 comparable 
with the size of absorbers as measured by their Doppler pa- 
rameter (10), (l + z)(d v )/l v ~ 3 - 10~ 2 [4/(1 ±z)] 1/2 . These 




0.001 0.010 

q, q/r 



0.100 



Figure 4. Top panel: the cumulative PDF, W t, s (q/T), together 
with fit (52) (long dashed line). Bottom panel: the correlation 
function of initial velocity field, 1 — together with the the- 

oretical fit (20) (dashed line) and the fit (60) (solid line) . Points 
and stars show the same functions found for the sample of 14 
QSOs and for the sample of 7430 separations. 



factors lead to the artificial blending of close absorbers, dis- 
tort their characteristics and decrease reliability of the esti- 
mates (59) for q r d < 10~ 2 . 

4-2.2 Correlation functions derived from DM column 
density of absorbers 

The same correlation function of the initial velocity field, 
£„((/), can be also found from estimates of the DM column 
density of absorbers, q. Here we use a more complex pro- 
cedure of determination of q from the observed z, b, & Nhi 
and poorly known Ti2 (57) what decreases its reliability. 
On the other hand, comparison of the functions found 
with two different approaches allows us to test the model of 
absorbers discussed in Sees. 3 and 4.1 . 

The observed cumulative PDF, W bs{q/r), and the re- 
constructed correlation function 1 — £„ (q) are plotted in Fig. 
4 for the samples of 19 QSOs with 6 270 and 7430 absorbers 
and 14 QSOs (3 674 absorbers). For these samples we have, 
respectively, 

W obs (< qM/W q (C fit ) w 1. ± 0.07, 0.9 ± 0.2 , 



i - £/«(«) 
i -£/«(?) 



0.95 ±0.15, 1.1 ±0.2, 
1.5q 



(60) 



p q = 0.8 • 10" 



^2.25g 4 ±g 2 ±pi'V' 6 

10~ 4 < q < 0.3, 3ft _1 kpc < l v q < 9.4ft _1 Mpc , 

where £/« = ((£fit,r) is given by (21). 

As is seen from Fig. 4, at q < 5 • 10~ 3 the derived 
correlation function (60) differs from the standard CDM - 
like one (20) . For these q both samples of absorbers are 
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quite representative with N aba w 1590 and N aba w 850 for 
samples of 19 and 14 QSOs, respectively. However, as was 
noted above, a more complex procedure of determination of 
decreases its reliability. 
At these scales the difference between the derived and 
reference correlation functions, (, v (q) and £cdm (20), could 
be mainly caused by the deficit of weaker absorbers (Sec. 3.9) 
because 1 — £„(g) oc W q 1 ^ 2 . As before, to test the impact of 
this factor we calculated the function (, v (q) for the sample 
of 7430 absorbers with lg Nhi < 15. For such sample the 
difference between the derived and reference functions, (q) 
and £,cDM{q), at q < 10~ 3 decreases and the correlation 
function is fitted by the expression 



!-£/«(«) = 



1.5<f 



p f = 4 • 10- 3 . (61) 



The functional forms of expressions (20), (60) and (61) are 
identical and they differ only by the values of fit parame- 
ters p q « 0.07po, Pf ~ 0.3po- This means that the differ- 
ence between (60) and £cdm (20) could be mainly re- 
lated to a possible incompleteness of the observed samples 
of weaker absorbers, to the limited precision of measure- 
ments of & & Nhi and to the limited precision of our model 
in describing such absorbers (see discussion in Sec. 3.7). 

The correlation function (61) is quite similar to the 
reference one and their difference is in the range of obser- 
vational errors. This fact demonstrates that probably the 
CDM like initial power spectrum can be traced at least 
down to q ~ 10~ 4 , l v q ~ 3/i _1 kpc. This means also that 
qo < 10~ 4 , mo > 25 and the effective mass of the DM parti- 
cles Mdm > 100 MeV. However, the function (61) is derived 
from inhomogeneous samples and therefore its reliability is 
in question. 

At larger scales, q > 5 • 10~ 3 , l v q > 0.15/i _1 Mpc, re- 
sults obtained from the analysis of both characteristics of 
absorbers are quite similar: 

1 - £v(qrd) 



1 — CcCAf(?rd) 



1.08 ±0.25, 1.06 ±0.17, 



1.0 ±0.2, 1.09 ±0.16. 



(62) 



(63) 



1 - ^CDM(q) 

These results confirm the CDM - like type of the initial 
power spectrum down to scales l v q > 0.15/i _1 Mpc what ex- 
tends conclusions of Croft et al. (2002), Viel et al. (2004b); 
McDonald et al. (2004) and Zaroubi et al. 2005. It also 
demonstrates the self consistency of the adopted model of 
absorbers. 

It is important, that similar results are found for both 
samples of 19 QSOs with 6 270 absorbers and of 14 QSOs 
with 3 674 absorbers. The samples used are compiled from 
spectra observed with different instruments and resolutions 
and the parameters of absorbers were found with different 
codes what increases their possible non homogeneity. Under 
these conditions, the stability of our results demonstrates 
their objectivity. 

Reconstruction of the initial power spectrum with the 
help of the relations (19) shows that in the range of errors 
the measured and CDM-like power spectra are quite simi- 
lar each other. At larger k, k/ko > 100, there is some excess 
of the power but estimates become unstable because of the 
limited range of measured q. Investigations of this impor- 



tant problem should be continued with more homogeneous 
sample of observed spectra. 



4.3 Statistical characteristics of the full sample of 
absorbers 

For the sample of 6 270 absorbers the redshift variations 
of the mean DM column density, (£), the mean fraction of 
DM component accumulated by absorbers, (f a bs), and the 
real and redshift sizes of absorbers along the line of sight, 

3/2 3/2 

{d a bsZ 4 }, and {d v z 4 ), are plotted in Fig. 5. The PDFs 
for the functions £ and d aba z^J 2 are plotted in Fig. 6. For ~ 
1 500 absorbers of this sample the parameter q was corrected 
as described in Sec. 3.7 (Eq. (49)). As was expected, the 
fraction of such absorbers (~ 25%) is close to the probability 
< cosf? < 0.25 to find absorbers oriented along the line of 
sight. 

The mean DM column density of absorbers, (q) or (£), 
is the most stable characteristic of the sample. In princi- 
ple, (£) does not change due to the formation and merging 
of absorbers and due to their transverse compression and/or 
expansion (DD04). It depends upon the ionization rate, Ti2, 
the amplitude of initial perturbations, To or <r 8 , and upon 
the shape of the correlation function of initial perturba- 
tions, £„, or the initial power spectrum, p(k). It also weakly 
depends upon the parameter bthr used to discriminate ab- 
sorbers formed by the adiabatic and shock compression. The 
differences between the expected and measured (f ) charac- 
terize, in fact, the completeness and representativity of the 
samples, the scatter of the function Ti2 and the influence of 
disregarded factors such as Q bg , Os & O x - 

For the sample under consideration the redshift vari- 
ations of the measured (£(«)) around the mean values are 
moderate, 



(C) « 0.82 ± 0.08, {qzl) w (1.8 ± 0.2) • 10" 



(64) 



At small and larger redshifts, z < 2 and z > 3.7, lim- 
ited statistic of absorbers decreases reliability of our esti- 
mates. The measured (C(z)) is close to the theoretically ex- 
pected value (25) what verifies the choice of the ionizing 
rate Ti2(z) in (57) and the amplitude of initial perturba- 
tions, ro = 0.19, cts ~ 0.9 . 

The PDF of the DM column density, P(() plotted in 
Fig. 6 is fitted with a scatter ~ 10% by the function 

P(x c ) = 1.2 exp(-a: c )erf {y/xfi/y/xc, x ( = C(«, z)/{Q , (65) 

which is very close to the theoretical relation (25) . This 
result verifies the self consistency of the physical model used 
here and the assumed Gaussianity of initial perturbations. 

As is well known, for the full sample the HI column 
density and Doppler parameter are weakly correlated, and 
for our sample their linear correlation coefficient is 



Rwi = \(bN HI ) - (b)(N HI )]/a b a HI « 0.16 , 



(66) 



where N^j = Nhi/z\. At the same time, the DM column 
density, £, is correlated with both the HI column density and 
Doppler parameter, and the linear correlation coefficients 
defined in the same manner as (66) are 

i? Cb «0.34, R ChI ^ 0.72. (67) 

The PDF of observed absorber separations, P{x r d), 
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Figure 5. Functions £(z)/(£}, (top panel), the matter frac- 
tion accumulated by absorbers, f ab s(z), (middle panels), and 

the real and redshift sizes of absorbers, xj (70), and x v = 

3/2 3/2 

d v (z)z 4 /(d v z 4 }, (bottom panel, solid and dashed lines). 




Figure 6. PDF for the DM column density, P(x(), (top panel), 
the absorber separation, P(x rc i), (middle panel), and for the ab- 
sorber thickness, P(x l {) (bottom panel) arc plotted together with 
fits (65), (68) and (70). 

plotted in Fig. 6 is fitted with a scatter ~ 13% by the func- 
tion 

P{x r d) = 3.56exp( — 1.6:r r d)erf 4 (VOiaVd)/v / aVd > (68) 

where x rd = ( r d(d sep , z)/{C, r d), (Crd) ~ 1-5 and Crd{d sep ,z) 
was introduced in (33). At small and large x r d this fit differs 
from the theoretically expected one for redshift space (33) 
but it is quite similar to the fit (30) for the PDF of separa- 



tions in Lagrangian space. Perhaps, this fact can be related 
to the unexpectedly moderate influence of the peculiar ve- 
locities of absorbers. However, it can be partly caused by 
peculiarities of our samples. 

Our model of absorbers (see Sec. 3) allows one to es- 
timate roughly the real size of baryonic distribution across 
the absorber, dabs, (41) . For the full sample we have 

(dibs) = {d abs zl /2 ) » (0.11 ± 0.01)fc _1 Mpc , (69) 

(d* v ) = (d v zl /2 ) w (0.12±0.01)/i _1 Mpc, 

and both sizes increase with time. The PDF of the real size 
of baryonic pancakes is roughly fitted by: 

P{x d ) w 0.6 exp(-0.8a> d ) + 0.4 cxp[-(0.73 - x d ) 2 , (70) 

x d = d abs /(d aba ) = d aba z 3 J 2 / (d aba z Z J 2 ) , 

with a scatter ~ 12%. 

For the expected mean transverse size of absorbers (27) 
and for model parameters (57) we have for the mean proper 
size 

(l v qtr)(l + z)- 1 « 0.45zJ 2 (l + 2)- 1 /i _1 Mpc , 

what is consistent with recent direct estimates by Becker, 
Sargent & Raugh (2004) at z ~ 3 - 3.5 

l v qtr ~ (0.15-0.2)/i _1 Mpc. 

The exponential PDF of the transverse sizes (27) and its 
strong redshift dependence explain large scatter of the sizes 
{l v qtr < lh~ Mpc) measured in many observations of pairs 
of QSOs (see, e.g., discussion in Becker, Sargent & Raugh 
2004). 

The mean measured fraction of matter accumulated by 
absorbers is 

(f ab s(z)) ^0.44 ±0.07. (71) 

In spite of the limited applicability of the one dimensional 
approach (31) and the limited precision of our model of ab- 
sorbers the measured fraction (71) is close to the theoreti- 
cal expectation of the Zel'dovich theory (32) what verifies 
the choice of the model characteristics (Sec. 4.1). The weak 
redshift variations of this function agrees well with the self 
similar evolution of absorbers. 

4.4 Adiabatically and shock compressed absorbers 

To investigate the complicated evolution of absorbers in 
more details, we compare subpopulations of adiabatically 
and shock compressed absorbers. These subpopulations were 
separated by comparison of the measured Doppler param- 
eter, b, with the background one (b bg ), (3). By definition, 
absorbers with b > b thr — 1.5(b bg ) belong to the subpopu- 
lation of shock compressed and strongly relaxed absorbers, 
while absorbers with b < b t hr = 1.5{b bg ) are considered as 
formed in a course of adiabatic or weak shock compression. 
This discrimination is not strict however and characteris- 
tics of subpopulations depend upon the sample used in the 
analysis and the parameters of the parameters of the back- 
ground (3, 57). This classification allows one to characterize 
reasonably well the observed absorbers and to trace their 
evolutionary history. For both subpopulations, redshift vari- 
ations of the mean characteristics are listed in Table 2 and 
some of the PDFs are plotted in Fig.7. 
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Table 2. Mean parameters of adiabatically and shock compressed 
absorbers 



adiabatic shock 



(In Ntt t /y 2 \ 

\ l 9 n hi 1 z 4/ 


1 Q o 4- n 07 


1 Q K _l_ n i 

ro.o in u.j- 


(b)km/s 


21 ±0.9 


48 ±5 


<fn) 


0.8 ±0.06 


0.2 ±0.05 


ifabs) 


0.33 ± 0.05 


0.11 ±0.08 


(0 


0.77 ±0.07 


1.1 ±0.2 


R( b 


0.39 


0.42 




0.75 


0.77 


(h) 


/3 7 / 2 


(2±0.3)z~ 2 ' 5 


(5 6 + 21n z 4 ) 




2.3 ±0.2 


(d o6i *2 /2 > h " lM P c 


0.09 ±0.01 


0.18 ±0.05 


<d t ,2f /2 ) /l ~ lM P C 


0.1 ±0.05 


0.22 ±0.02 



Z4 = (1 + z)/4, f n and f a b a are the fraction of absorbers and of 
matter in absorbers, i?<* &cRq h i are the linear correlation coef- 
ficients of hydrogen and DM column densities and Doppler pa- 
rameter defined in the same manner as (66), <5(, and Sb are the 
overdensity above the mean density and relative entropy of com- 
pressed baryonic component, d a bs and d v are the real and rcdshift 
sizes of absorbers along the line of sight. 

For the sample under investigation ~ 80% of absorbers 
are compressed adiabatically and they accumulate ~ 80% of 
the compressed matter. These fractions, f n &zf a bs, weakly 
vary with redshift. The PDFs, P(x c ),x c = C/(C>, plotted 
in Fig. 7 for both subpopulations are quite similar to each 
other and to the PDFs (34) and P(x ( ) (65) plotted in Fig. 
6 for the full sample. It demonstrates that in wide range 
of redshifts absorbers could be formed by both processes 
while the cutoff at f = 0.3(f) in the PDF for shock com- 
pressed absorbers is imposed by the method used for the 
discrimination of absorbers. For both samples, the correla- 
tion coefficients R( b and R( H i defined in the same manner 
as (66) are similar to each other. These similarities verify 
the generic nature of absorbers and indicate that they can 
be successfully combined into one sample, what is consistent 
with expectations of the Zel'dovich theory. For the subpop- 
ulation of shock compressed and strongly relaxed absorbers 
(C) ~ {(rnrg) (34), as expected for merged absorbers, what 
confirms the importance of merging in the process of ab- 
sorbers' evolution. 

For both subpopulations (b/b in f) < 1, what can be 
partly related to evaporation of high velocity particles in 
the course of relaxation of the compressed matter. For both 
subpopulations, the redshift size of absorbers along the line 
of sight, d v , is larger than the real size, d a b s (Table 2). The 
redshift dependence d a bs oc z~ 3 /' 2 is consistent with the evo- 
lution of their redshift size. The PDFs of the real size of 
absorbers, P(xd), is plotted in Fig. 7. 

For adiabatically compressed absorbers the entropy is 
the same as for the background while the overdensity, S b = 
/3 3 , is determined by b and bb g . In particular, for cold ab- 
sorbers with b < btg we have Sb < 1. 

For subpopulation of strongly relaxed and shock com- 
pressed absorbers both the entropy and the overdensity de- 
pend upon the complex evolutionary history of absorbers 
and cannot be described by a simple theoretical model. For 
such absorbers growth of the overdensity with time can be 
explained by the action of several factors such as successive 
merging and contribution of long lived absorbers formed at 
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Figure 7. The PDFs for the DM column density, P(x^), overden- 
sity, P(xs), and the real size across the absorber, P(x£), (three 
top panels) arc plotted for shock compressed (solid lines) and adi- 
abatically compressed (dashed lines) absorbers together with fits 
(65) and (72) . For shock compressed absorbers, the PDFs for the 
baryonic entropy, P(xf), is plotted in bottom panel with fit (73). 

redshifts Zf larger than the observed one, z a b s < Zf . For 
such absorbers the measured PDF P(x ) is plotted in Fig. 
7 and it can be fitted by a superposition of two exponential 
functions: 

P(x s ) « 3 exp(-3.1x s ) + 0.13 exp(-0.2^ 5 ) , (72) 
x s = 8b(l + zf- 5 /{8 b {l + z) 2 - 5 ) . 

The entropy of strongly relaxed and shock compressed 
baryons (50) increases with time and its PDF plotted in 
Fig. 7 is well fitted by the Gauss function 

P(x f ) w exp[-(x/ - 0.95) 2 /0.35] , (73) 

x f = \n[F s (z)zl]/(\n[F s (z)zt]) . 

Such PDF naturally arises when the entropy is generated by 
the action of many random factors such as the shock waves 
accompanying the successive merging of absorbers. 

4.5 Absorbers in rapidly and moderately 
expanded regions 

As was discussed in Sec. 3.5, we expect that majority of 
absorbers with small Doppler parameters, b < b rap ^23.5 
km/s, and especially with b < (bb g ) could be formed within 
rapidly expanded regions and, so, they can characterize some 
properties of these regions. Some characteristics of absorbers 
with b < b rap are listed in Table 3 in comparison with the 
same characteristics of adiabatically compressed absorbers 
with b t hr > b > b rap situated within moderately expanded 
regions. 

The rapidly expanded regions accumulate (/„) ~ 50% 
of adiabatically compressed absorbers and (frap) ~ 30% of 
adiabatically compressed matter. These values practically 
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Table 3. Parameters of absorbers in rapidly and moderately ex- 
panded regions 





rapidly 


moderately 


(IgNmfzl) 


13. ±0.1 


13.4 ±0.1 


(b)km/s 


15.5 ±0.6 


27 ± 1.2 


<fn) 


0.49 ±0.05 


0.51 ±0.05 


{fabs ) 


(0.11 ± 0.05) 


(0.22 ±0.05) 


(0 


0.46 ±0.04 


1 ± 0.1 




0.8 


0.07 


Rq h i 


0.4 


0.76 


{D re gzl)h- 1 Mpc 


(1.6 ±0.2) 


(3.2 ±0.4) 


(d^z^h^Mpc 


(0.08 ± 0.01) 


(0.1 ±0.01) 


(d v zl /2 )h- 1 Mpc 


(0.08 ± 0.003) 


(0.12 ± 0.005) 



24 = (1 + z)/4, (f„) and (f a bs) are the fraction absorbers and 
of matter in absorbers, Rr, & Rr , are the linear correlation 

4 

coefficients of hydrogen and DM column densities and Doppler 
parameter defined in the same manner as (66), (D reg ) is the mean 
size of rapidly and slowly expanded regions, (d sep ) and (d v ) are 
the mean absorber separation and their mean redshift size along 
the line of sight, 



4.6 Absorbers and properties of the background 

Some characteristics of absorbers could be used to estimate 
the redshift variations of the mean properties of homoge- 
neously distributed hydrogen (see, e.g., Hui & Gnedin 1997; 
Schaye et al. 1999, 2000; McDonald et al. 2001). However, 
such estimates are inevitably approximate and their signifi- 
cant scatter is caused by the action of many random factors 
discussed above. To obtain more stable results Schaye et al. 
(1999, 2000) consider a cutoff at small b in the distribution 
of b(Niji)- However, such absorbers are probably formed 
within rapidly expanded regions and for them both back- 
ground properties and expansion rate vary randomly from 
absorber to absorber. 

For the subpopulation of moderately expanded ab- 
sorbers, we can combine Eqs. (41), (42), (36) and (46) and, 
in principle, connect the background temperature with the 
Doppler parameter of absorbers, b, and their hydrogen col- 
umn density, Nhi- However, reliability and significance of 
such estimates are in question. Some restrictions on the in- 
tensity of UV background were discussed in Sec. 4.1 . 



do not depend on the redshift. The Doppler parameter, b, 
and the DM column density, ( rap (q, z), also weakly depend 
upon redshift. For this subpopulation we have 

(Crap) ~ 0.5(C), (frap) ~ 0.25{f abs ) . (74) 

As was noted in Sec. 2.2, the mean size of rapidly ex- 
panded regions, (D rap ), as well as the mean separation of 
absorbers within these regions, {d sep ), increase with time in 
the same manner as the mean separation of absorbers for 
the full sample (10). At zs, ~ 1 the typical mass associated 
with rapidly expanded regions, M rap (14), is in the range of 
galactic masses and it increases with time oc (1 + z)~ 3 . This 
result is consistent with the expected symmetry of positive 
and negative initial density perturbations what leads to for- 
mation of both galaxies and rapidly expanded regions. These 
results are consistent with theoretically expected ones (37), 
what confirms the interpretation of the complex shape of 
PDF Pj, (6) and subpopulation of weak absorbers proposed 
in Sec. 3.5 . 

For the subpopulation of moderately expanded ab- 
sorbers variations of the Doppler parameter, b, axe small, 
R(b <SC 1 and the DM column density, C, depends mainly 
upon the hydrogen column density, Nhi- In contrast, for the 
rapidly expanded absorbers the influence of both b & Nhi 
are equally important. In spite of this difference, for both 
subpopulations the PDFs P(C) are quite similar to each 
other and to the PDF (65) obtained for the full sample. This 
fact indicates that the interaction of small and large scale 
perturbations changes (C) more strongly but only weakly 
influences the shape of the PDFs P(C)- 

For both subpopulations, the mean proper sizes of ab- 
sorbers, (d a bs), are similar (Table 3) but their PDFs are 
quite different. Thus, for rapidly expanded regions the PDF 
P{xd) is step-like and it is responsible for the bump at 
datsz 3 / 2 < {d abs zl /2 } in the PDF P(x d ) plotted in Fig. 6. 
This distribution differs from the distribution of Doppler pa- 
rameter, P(xb), what suggests a complex internal structure 
of such absorbers. 



5 SUMMARY AND DISCUSSION. 

In this paper we continue the analysis initiated in Paper I 
and Paper II that is based on the statistical description of 
Zel'dovich pancakes (DD99, DD04). This approach allows 
one to connect the observed characteristics of absorbers with 
fundamental properties of the initial perturbations without 
any smoothing or filtering procedures, to reveal and to illus- 
trate the main tendencies of structure evolution. It demon- 
strates also the generic origin of absorbers and the Large 
Scale Structure observed in the spatial distribution of galax- 
ies at small redshifts. 

We investigate the more representative sample of 
~ 6 000 absorbers what allows us to improve the physical 
model of absorbers introduced in Paper I and Paper II and 
to obtain reasonable description of physical characteristics 
of absorbers. The progress achieved demonstrates again the 
key role of the representativity of the observed samples for 
the construction of the physical model of absorbers and re- 
veals a close connection between conclusions and the ob- 
servational database. Further progress can be achieved with 
richer and more refined sample of observed absorbers. 

5.1 Main results 

Main results of our analysis can be summarized as follows: 

(i) For suitable parameters of the model (Sec. 4.1), the 
basic observed properties of absorbers and their evolution 
are quite successfully described by the statistical model of 
DM confined structure elements (Zel'dovich pancakes) with 
various evolutionary histories. Comparison of independent 
estimates of the DM characteristics of pancakes confirms 
the self consistency of the physical model. This model is in a 
good agreement with measured properties of metal systems 
(see, e.g. Carswell, Schaye & Kim 2002; Telferet al. 2002; 
Bergeron & Herbert-Fort 2005). 

(ii) The PDFs of the DM column density and the dis- 
tances between neighboring absorbers are found to be consis- 
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tent with the Gaussian initial perturbations with the CDM- 
like initial power spectrum. 

(iii) For the observed range of redshifts the evolution of 
absorbers is close to self-similar one. This implies that it 
leads to slow variations of mean absorber characteristics 
with redshift and retains their PDFs. 

(iv) We estimate the shape of the correlation function of 
the initial velocity field what in turn allows us to estimate 
the shape of the initial power spectrum. At scales > 0.15/i -1 
Mpc both derived correlation functions, (59) and (60), repro- 
duce the CDM-like one. This means that at such scales the 
power spectrum of initial perturbations is close to the stan- 
dard one (16) . At smaller scales we see some differences be- 
tween the derived and CDM-like correlation functions which 
depend upon the sample used in the analysis. 

(v) Analysis of variations of the Doppler parameter, 6, 
along the line of sight demonstrates existence of rapidly 
expanded regions which can be considered as examples of 
strong negative density perturbations of galactic mass scale. 

(vi) Our analysis shows that in the observed range of 
redshifts we can expect slow variations of the intensity of 
UV background radiation and of the ionization rate, Ti2. 
Our results are close to the estimates of the UV back- 
ground in Haardt&Madau (1996), Scott et al. (2002), and 
Demianski & Doroshkevich (2004b). 

5.2 Test of the model of absorbers 

The physical model of absorbers introduced in Sec. 3 links 
the measured z, b and Nhi with other physical character- 
istics of both gaseous and DM components forming the ob- 
served absorbers. It is important that this ID model pro- 
vides us with the self consistent statistical description of the 
Ly-a forest although some parameters of pancakes remain 
unknown. Action of these parameters as well as uncertain- 
ties in the available estimates of the background temper- 
ature and the UV background radiation lead to moderate 
random scatter of the derived characteristics of absorbers. 
Fortunately, actions of these factors partly compensate each 
other, what allows us to obtain reasonable statistical de- 
scription for majority of absorbers. The self consistency of 
this approach is confirmed by similarity of the functions 
(,v{q r d) (59) and £ v (q) (60) and by estimates of the matter 
fraction, f a bs, accumulated by absorbers (71). These func- 
tions are related to independent characteristics of absorbers 
obtained from measurements of their separation and their 
DM column density. 

For richer absorbers, both the hydrostatic equilibrium 
of compressed matter along the shorter axis of pancakes and 
the close link between the gas temperature and the Doppler 
parameter, b, are confirmed by comparison of characteris- 
tics of the HI and metal systems (see, e.g. Carswell, Schaye 
& Kim 2002; Telferet al. 2002; Simcoe, Sargent & Rauch 
2002, 2004; Boksenberg, Sargent and Rauch 2003; Manning 
2002, 2003 a,b; Bergeron & Herbert-Fort 2005). In particlar, 
for 191 high resolution metal systems presented in Boksen- 
berg, Sargent and Rough (2003) differences between the gas 
temperatures measured by the Doppler parameters of HI 
and CIV do not exceed ~ 25%, what is comparable to the 
precision of measurements. Comparison of the Doppler pa- 
rameters measured for HI, CIV and OVI (Carswell, Schaye 
& Kim 2002) verifies also their similarity and shows that as 



a rule the macroscopic (turbulent) velocities are subsonic. 
These observational results strongly support the domina- 
tion of long-lived gravitationally bound and partly relaxed 
absorbers composed of both DM and baryonic components. 

Numerical simulations show that the line width depends 
upon the thermal broadening, the differential Hubble flow 
and peculiar velocities and the relative influence of these 
factors varies from absorber to absorber (see, e.g., Theuns, 
Schaye & Haehnelt 2000; Schaye 2001). The Hubble flow is 
more essential for weaker absorbers and can artificially in- 
crease their Doppler parameter. For majority of absorbers in 
the considered ID model the possible contribution of Hubble 
flow is naturally linked with the compression or expansion of 
pancakes in the transverse directions and depends upon the 
(unknown) relative orientation of absorber and the line of 
sight. The available observational data do not allow to dis- 
criminate between the thermal and macroscopic broadening 
of the lines what increases the random scatter of our results. 
To perform such discrimination a more detailed description 
of the observed line profiles is required. 

One of the important problem facing the high resolution 
numerical simulations is the development of the methods for 
the more detailed reconstruction of the physical properties 
and revealing of links between the DM and baryonic compo- 
nents of observed absorbers. In particular, this includes the 
discrimination of the thermal and macroscopic broadening 
of lines, explanation of the surprisingly weak redshift de- 
pendence of the mean Doppler parameter, detection of the 
complex internal structure of absorbers as indicated by the 
observations of metal systems and so on. 

However, now technical limitations restrict facilities of 
simulations. Thus, the small box size used eliminates the 
large scale part of the power spectrum and decreases the rep- 
resentativity of simulated sample of absorbers. As was dis- 
cussed in Paper II and in Manning (2003 a,b), these factors 
eliminate the interaction of large and small scale perturba- 
tions and distort characteristics of the simulated absorbers. 
Simulations reproduce the observed transmitted flux and its 
main features and now they are used mainly for the sur- 
prisingly stable reconstruction of the initial power spectrum 
from the flux characteristics (see, e.g., Seljak et al. 2004; Mc- 
Donald et al. 2004, Viel et al 2004a, b). However, the analysis 
of Meiksin, Bryan and Machacek (2001) shows that simula- 
tions have problems with reproduction of the observed PDFs 
for the column density of neutral hydrogen, NhI, and the 
Doppler parameter, b, and their self similar redshift evolu- 
tion. 

More detailed criticism of the "Fluctuating Gunn- 
Petersen approximation" and the simulations of the forest 
can be found in Manning (2002, 2003 a,b), Paper II and 
references cited in these papers. 

5.3 Properties of absorbers 

Analysis of the mean absorbers characteristics performed 
in Sec. 4 shows that the sample of observed absorbers is 
composed of pancakes with various evolutionary histories. 
We discuss five main factors that determine evolution of 
absorbers after their formation. They are: the transverse 
expansion and compression of pancakes, the disruption of 
structure elements into a system of high density clouds, the 
merging of absorbers and the radiative heating and cooling 
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of compressed gas. The first two factors change the overden- 
sity of DM and gas but do not change the gas entropy. Next 
two factors change both the gas entropy and overdensity but 
do not change the DM characteristics. 

The sample of observed forest can be naturally divided 
into subsamples of adiabatically and shock compressed ab- 
sorbers formed by merging. Moreover, about half of adia- 
batically compressed absorbers are formed within rapidly 
expanded regions where the background temperature is less 
than mean one. So, the temperature of absorbers formed 
within such regions can be also less than the mean tem- 
perature of the background (3). These results illustrate the 
influence of some of the factors mentioned above. However, 
the slow variation of the mean characteristics of absorbers 
and their PDFs with redshift confirms that we observe the 
self-similar period of absorbers evolution when the action of 
these factors is balanced and regular variations of the UV 
background does not distort this balance. 

For shock compressed absorbers, introduction of the 
DM column density, q and £ , and entropy, S m ~ St, allows to 
discriminate between the systematic and random variations 
of their properties. The former ones are naturally related to 
the progressive growth with time of the DM column density 
of absorbers, q(z)&z£(z), and they can be described theo- 
retically. On the other hand the action of random factors 
cannot be satisfactorily described by any theoretical model. 
However, in the framework of our approach, the joint action 
of all random factors is summarized by one random func- 
tion, St, directly expressed through the observed parame- 
ters (50). These results alleviate the problem of description 
of absorbers and, perhaps, the modelling of the Ly-a forest 
based on the simulated DM distribution (Viel et al. 2002) . 

For adiabatically compressed absorbers, the spatial dis- 
tributions, entropy and overdensity of baryonic and DM 
components are different. Unfortunately at present these 
characteristics of DM component cannot be determined from 
observations with a reasonable accuracy. For such absorbers 
the baryonic entropy is identical to the background one given 
by (4) while the PDFs and the random scatter of observed 
characteristics are determined mainly by random variations 
of the expansion rate and the background density and tem- 
perature. For this subpopulation, the process of formation 
and evolution of absorbers should be investigated more thor- 
oughly. 

5.4 Characteristics of the initial power spectrum 

The initial power spectrum of density perturbations is cre- 
ated at the period of inflation and its observed determination 
is very important for investigations of the early Universe. 
The amplitude and the shape of large scale initial power 
spectrum are approximately established by investigations of 
relic radiation (see. e.g, Spergel et al. 2003, 2006) and the 
structure of the Universe at z < 1 detected in large redshift 
surveys such as the SDSS (Dodelson et al. 2002; Tegmark 
et al. 2004) and 2dF (Percival et al. 2001). The shape of 
the initial power spectrum at small scale can be tested at 
high redshifts where it is not so strongly distorted by non- 
linear evolution (see, e.g., Croft et al. 2002; Tegmark et al. 
2002; McDonald et al. 2004; Seljak et al. 2004; Zaroubi et 
al. 2005). 

Here we retrieve the correlation function of initial ve- 



locity field, f„, from direct measurements of the PDFs of 
fundamental characteristics of absorbers such as their sepa- 
ration, d sep , and the DM column density, q . Both estimates 
are derived in the same way and result in the same shape of 
the correlation function at larger scales, q > 5 • 10~ 3 ,Z„g > 
0.15ft -1 Mpc . At these scales the measured correlation func- 
tions coincide with the CDM-like one (20) what confirms 
conclusions of Croft et al. (2002), Viel et al. (2004b); Mc- 
Donald et al. (2004) and Zaroubi et al. (2005) obtained at 
scales > lh^ 1 Mpc. 

At smaller scales the results obtained with analysis of 
the absorbers separation, d sep , and the DM column density, 
q, are different and demonstrate some excess of power at 
scales 150 kpc > l v q > 3 kpc. Parameters of these functions 
and the excess depend upon the sample of absorbers and for 
the extended samples the derived correlation functions be- 
come quite similar to the CDM-like one. However, reliability 
of this result is in question due to a probable incompleteness 
of the extended samples. 

The interpretation of these distortions is not unique be- 
cause of very limited available information. As was shown in 
Sec. 4.2, they are sensitive to the deficit of weaker absorbers 
and small separations of absorbers in the sample under con- 
sideration. Therefore these distortions can be enhanced by 
the probable incompleteness of the observed sample created 
by the finite resolution of observations, blending of lines and 
approximate character of our analysis. If this explanation 
is correct than these factors restrict the presently available 
range of investigations to l v q > 100 kpc. Further progress 
can be achieved with more refined observations of absorp- 
tion spectra of QSOs and with more refined identification of 
absorbers in the observed spectra. 

In turn, these divergences can be related to special fea- 
tures in the initial power spectrum at small scales. Recent 
WMAP measurements indicate that adiabatic Gaussian per- 
turbations dominate on large scale (Peiris et al. 2003; Ko- 
matsu et al. 2003). However, these results do not preclude 
deviations from the standard CDM-like power spectrum at 
small scales. In particular, such deviations appear in models 
of the one field inflation with a complicated inflation poten- 
tial (see, e.g., Ivanov, Naselsky & Novikov 1994) or multiple 
fields inflation (see, e.g., Polarski & Starobinsky 1995; Turok 
1996). Both models generate adiabatic or isocurvature de- 
viations from the simple CDM-like power spectrum. More 
detailed discussion of such models can be found, for exam- 
ple, in Peiris et al. (2003). 

5.5 Absorbers as elements of the Large Scale 
Structure of the Universe 

At redshifts z > 1.7 the Large Scale Structure is observed 
mainly as systems of absorbers in spectra of distant QSOs. 
Numerical simulations show that even at such redshifts we 
can see also high density filaments and clumps formed by 
"galaxies" and some of them are actually observed in spec- 
tra of QSOs as metal systems, Lyman damped and Lyman 
limit systems. However, available observational data do not 
yet allow one to characterize statistically properties of such 
structure elements (see, e.g., Boksenberg, Sargent & Rauch 
2003). 

At small redshifts, the Large Scale Structure is observed 
as a spatial distribution of both galaxies and neutral hydro- 
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gen. The investigation of galaxy distribution in the SDSS 
DR1 (Doroshkevich, Tucker, Allam & Way 2004a) results in 
estimates of typical parameters of galaxy walls as 

(q) « 0.4, (b) « 320km/ s, {d sep } « 60/i _1 Mpc . (75) 

With these data the expected column density of neutral hy- 
drogen within the typical wall (38) is Nhi ~ 10 11 cmT 2 and 
even so spectacular object as the 'Greet Wall' does not man- 
ifest itself through absorbers. 

Our results indicate the generic link of absorbers and 
DM Zel'dovich pancakes and demonstrate that the embryos 
of walls could also be seen already at z ~ 3. Indeed, for 
basic parameters of subpopulation of 1 370 shock compressed 
absorbers with b > 30 km/s we have 

(d sep ) « (50 ± 11)(1 + zy 2 h- 1 Mpc , (76) 
(q) « (0.4±0.07)(l + z)" 2 , 

what is quite similar at z~0 to that given in (75). This fact 
indicates that, in principle, such absorbers can be consid- 
ered as embryos of wall-like elements of the Large Scale 
Structure of the Universe. Of course, such identification of 
walls observed in the galaxy distribution with elements of 
Ly-a forest is quite arbitrary and ignores the actual com- 
plex evolution of the LSS elements. However, it confirms 
generic character of the LSS evolution from richer absorbers 
to galaxy walls. The problem deserves further investigation 
first of all with more representative numerical simulations. 

For the first time poor absorbers at small redshifts were 
observed by Morris et al. (1991, 1993) and ~ 1000 of such 
absorbers were found by Bahcall et al. (1993, 1996) and Jan- 
nusi et al. (1998). Some of these absorbers are identified with 
halos of galaxies (see, e.g., Lanzetta et al. 1995; Le Brune, 
Bergeron & Boisse 1996) or galaxy filaments (Penton, Shull 
& Stock 2002) but others are situated far from any galaxies. 
These observations demonstrate that the space between the 
LSS elements - so called 'voids' - is not empty and contains 
essential fraction of baryonic and DM components of the 
matter. 

More detailed characteristics of absorbers at z <C 1 are 
given in Penton, Shull & Stock (2000, 2002); McLin et al. 
(2002) where the main absorber properties are found to be 
similar to those observed at high redshifts. For 79 absorbers 
with 12 < Ig Nhi < 15, 11km/ s < b < 80km/s listed in 
these papers the mean absorber separation is 

(d sep ) ~ (10±3)ft- 1 Mpc. (77) 

Despite the strong difference of many conditions at z < 1 
and z > 1.5, these observed characteristics of absorbers are 
quite similar to expected ones (10) extrapolated to z — 0. 

5.6 Observed and expected evolution of the Large 
Scale Structure 

Comparison of the expected and derived from observations 
characteristics of absorbers demonstrates that at the ob- 
served range of redshifts, 1.7 < z < 4.5 we see the self - 
similar period of structure evolution. During this period the 
main factors determining the evolution of absorbers such as 
the pancake expansion, creation and merging, are balanced 
what leads to relatively slow evolution of the mean proper- 
ties of absorbers, such as d S e P and the DM column density, 



q . This slow evolution is supported by slow regular vari- 
ations of the UV background radiation and the ionization 
rate, Ti2 . 

However, at small redshifts, z < 0.5 the growth of per- 
turbations and merging of absorbers becomes decelerated 
due to the influence of the A-term while expansion and dis- 
ruption of absorbers remains important. This means that at 
such redshifts the quiet evolution of absorbers is distorted 
and we can expect a progressive decrease of linear density 
of observed absorbers with the hydrogen column density 
Nhi > 10 12 cm~ 2 . The variations of the population of ob- 
served absorbers are also modulated by the poorly known 
variations of the UV background. 

At larger redshifts evolution of DM pancakes is mainly 
driven by the shape of the initial power spectrum. For the 
standard CDM-like correlation function (20) we get that the 
self-similar evolution takes place at redshifts 

z < z thr ~ v/0.375/go w 19^/10-3/90 , (78) 

and at z > z t hr pancakes with q ~ qo are more abundant. 

However, the observational test of these expectations is 
quite problematic because the observed characteristics of ab- 
sorbers depend also upon evolution of the background tem- 
perature and UV radiation. 



5.7 Reheating of the Universe 

Recent observations of high redshift quasars with z > 5 
(Djorgovski et al. 2001; Becker et al. 2001; Pentericci et 
al. 2002; Fan et al. 2002, 2003, 2004) provide clear evi- 
dence in favor of the reionization of the Universe at redshifts 
z ~ 6 when the volume averaged fraction of neutral hydro- 
gen is found to be fn > 10 -3 and the photo ionization rate 
T 7 ~ (0.02 - 0.08) • lO^s" 1 . These results are consistent 
with those expected at the end of the reionization epoch 
which probably takes place at z ~ 6. 

These results can be compared with expectations of 
the Zel'dovich approximation (DD04) . The potential of this 
approach is limited since it cannot describe the nonlinear 
stages of structure formation and, so, it cannot substitute 
the high resolution numerical simulations. However, it de- 
scribes quite well many observed and simulated statistical 
characteristics of the structure such as the redshift distribu- 
tion of absorbers and evolution of their DM column density. 
This approach does not depend on the box size, number of 
points and other limitations of numerical simulations (see 
discussion in Paper II) and it successfully augments them. 

This approach shows (DD04) that at z ~ 6 only ~ 3.5% 
of the matter is condensed within the high density clouds 
which can be associated with luminous objects. This value 
can increase up to ~ 5 - 6% with more accurate description 
of the clouds collapse. The same approach also allows one 
to estimate the mass function of structure elements (DD04) 
at different redshifts. At z ~ 6, the mean DM mass of the 
clouds is expected to be (M c i) ~ 10 10 Mq and majority of 
clouds have masses between 10 _3 (M ci ) and 10 {M d ). The 
formation of low mass clouds with M cl < 10 6 M Q is sup- 
pressed due to strong correlation of the initial density and 
velocity fields at scales < l p ~ 0.03ft _1 (g /10 _3 ) Mpc (17). 
However, the numerous low mass satellites of large cen- 
tral galaxies can be formed in the course of disruption of 
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massive collapsed clouds at the stage of their compression 
into thin pancake-like objects (Doroshkevich 1980; Vishniac 
1983). The minimal mass of such satellites was estimated in 
Barkana, Haiman & Ostriker (2001). 

This means that the investigation of absorbers observed 
at high redshifts should be supplemented by the study of 
properties of dwarf isolated galaxies and discrimination be- 
tween such galaxies and dwarf satellites of more massive 
galaxies. Such observations seem to be a perspective way to 
discriminate between models with one and several types of 
DM particles. 
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